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bath coupled by long-range hyperfine-mediated interactions 
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We investigate decoherence due to pure dephasing of a localized spin qubit interacting with 
a nuclear spin bath. Although in the limit of a very large magnetic field the only decoherence 
mechanism is spectral diffusion due to dipolar flip-flops of nuclear spins, with decreasing field the 
hyperfine-mediated interactions between the nuclear spins become important. We take advantage 
of their long-range nature, and resum the leading terms in an 1/N expansion of the decoherence 
time-evolution function (N, being the number of nuclear spins interacting appreciably with the 
electron spin, is large). For the case of the thermal uncorrelated bath we show that our theory 
is applicable down to low magnetic fields (~ 10 mT for a large dot with N = 10 6 ) allowing for 
comparison with recent experiments in GaAs quantum dot spin qubits. Within this approach we 
calculate the free induction decay and spin echo decoherence in GaAs and InGaAs as a function 
of the number of the nuclei in the bath (i.e. the quantum dot size) and the magnetic field. Our 
theory for free induction decay in a narrowed nuclear bath is shown to agree with the exact solution 
for decoherence due to hyperfine-mediated interaction which can be obtained when all the nuclei- 
electron coupling constants are identical. For the spin echo evolution we show that the dominant 
decoherence process at low fields is due to interactions between nuclei having significantly different 
Zeeman energies (i.e. nuclei of As and two isotopes of Ga in GaAs), and we compare our results with 
recent measurements of spin echo signal of a single spin confined in a GaAs quantum dot. For the 
same set of parameters we perform calculations of decoherence under various dynamical decoupling 
pulse sequences, and predict the effect of these sequences in low B regime in GaAs. 



I. INTRODUCTION 



The quantum dynamics of a localized central spin cou- 
pled to an environment of other fluctuating spins (spin- 
bath) is probably the longest studied model of quan- 
tum decoherence, outgrowing from research on dephasing 
and relaxation of ensembles of spins in Nuclear Magnetic 
Resonance!^! I n recent years this field of research has 
been reinvigorated by experimental progress in creation, 
initialization, readout, and coh erent control of electron 
spin qubits in semiconductors.^* 4 * 5 * 6 * 7 * 8 * 9 * 10 * 11 * 12 * 13 * 14 * 15 * 
Understanding the decoherence caused by the nuclear 
spin-bath is crucial for development of these qubits, since 
at low temperatures and moderate magnetic fields all 
the other sources of electron spin relaxation and de- 
phasing (e.g. spin- orbit scat tering with phonons) are 
strongly suppressed! 16 * 17 * 1 ^ 9 ^ The remaining relevant in- 
teraction is the hyperfine (hf ) interaction of an electron 
spin with the nuclear spins of the host material. In 
quantum computing architectures based on III-V ma- 
terials such as GaAs and InAs the presence of the nu- 
clear bath is unavoidable, since all the isotopes of Ga, 
As, and In have non-zero nuclear spinsJ^Sl The problem 
at hand is hard, since the bath is interacting: the nuclei 
are effectively interacting due to their mutual hf cou- 
pling to the "central spin" of the electron, and further- 
more they are also coupled by "intrinsic" intra-bath dipo- 



lar interactions! 20 * 21 ^ The resulting dynamics of the elec- 
tron is of highly nontrivial non-Markovian character,^ 
since time-scales of the bath dynamics are orders of 
magnitude longer than the typical electron spin preces- 
sion and decoherence time-scales. A simple per turba tion 
theory in hf interaction quickly breaks downJ22E2 anc j 
one has to use more global approaches (e.g. involving 
infinite-order resummations of the perturbation theory) 
in order to obtain physical results. However, significant 
theoretical progres s has been made with a combination 
of analytical Vv-orl J 21 ' 22 ' 23 ' 24 ' 25 ' 26 ' 27 ' 28 ' 29 ' 30 ' 31 ' 32 ' 33 ' 34 ' 351361 
and exact numerics for small systems^ 20 * 37 * 38 * 39 * 40 * 41 * 42 1 
A successful co mpar ison between the theory^J 27 * 28 * 30 * 4 ^ 
and experiment! 44 * 45 ! has been made for Si:P spin qubits 
(based on the spin of electron bound to the phosphorus 
donor). These calculations have shown that in experi- 
ments on Si:P system the decoherence is due to the dipo- 
lar interactions between the nuclei - the so-called spectral 
diffusion (SD) process, in which the hf interaction serves 
only as a transmitter of the dipolar-induced nuclear dy- 
namics to the electron. The analogou s calculations for 
GaAs dots were only 

performed 25]28129| for 

magnetic fields 

much higher than B«0.1T used in experiments on de- 
coherence of a single spin in a GaAs dotj^ and they 
gave time-scales two orders of magnitude longer than the 
observed 1 ^ decay in w 0.3 fis. This clearly showed that 
decoherence at low field in GaAs is not due to spectral 
diffusion. At fields currently experimentally relevant for 
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III-V based spin qubits one clearly has to consider the nu- 
clear dynamics caused by the hf interaction, specifically 
by the nuclei-nucle i interac tions due to virtual electron- 
nuclear spin flipsP^- 6 * 30 * 33 ! In this paper we present the 
detailed derivation of the theory^ of decoherence due to 
these hf-mediated interactions which is applicable down 
to ~ 10 (100) mT fields in quantum dots containing ~ 10 6 
(10 4 ) nuclei. 
The hf interaction is given by 

H M = £ A*S • J; = £ MS' J' + \{S + J^ + S-J?)} , 

i i 

_ (i) 

where A; is the hf coupling of the i-ih nuclear spin to 
the electron spin. This interaction causes both the en- 
ergy relaxation of the qubit (decay of the spin component 
parallel to the magnetic field) and its dephasing (decay of 
components of the spin transverse to the magnetic field) . 
The energy re laxation is strongly suppressed by a finite 
magnetic field j22123J24| j ue ^ l ar g e mismatch of elec- 
tronic and the nuclear Zccman energies the hf interaction 
is not efficient at relaxing the z component of the spin at 
fields larger than ~ 10 mT in GaAs quantum dots. More 
precisely, the longitudinal polarization loss was showrP^l 
to be of the order of S 2 , with the small parameter 5 de- 
fined as 



where A = ^ ^% is the total hf interaction energy, f2 is 
the electron spin Zccman splitting, and N is the number 
of nuclei interacting appreciably with the electron. In 
quantum dots the value of N ranges from 10 4 in sm allest 
self-assembled dotsPESl to 10 6 in lateral gated dots, 1121121 
and typical values of A and f2 are given in Table [T] 

The dephasing of the qubit, on the other hand, while 
being detrimental for quantum computation applications, 
is much harder to suppress. Theoretically, the simplest 
possible way to eliminate dephasing is to spin-polarize 
the nuclear bath, but the > 90% polarization ne cessary to 
achieve a meaningful increase in cohe rence ti mcP ^ 31 * 39 ! is 
impossible to achieve experimentally. ^ 48 * 49 ! Here we fo- 
cus on electron spin dephasing due to hf interaction with 
an unpolarized nuclear bath. Specifically, we consider 
a thermal uncorrelated state of the nuclei, which due to 
the smallness of the intrinsic nuclear energy scale (nuclear 
Zeeman splitting is ~ mK at B fields of the order of a 
Tesla) means that we will use the nuclear density matrix 
proportional to unity. Let us note that specially designed 
experimental protocols have been shown to drive the nu- 
clear bath into a different state, characterized by a small 
net polarization and, more importantly, a "narrowed" 
distribution of nuclear polarizations 

One operational definition of the spin dephasing time 
is the characteristic time of Free Induction Decay (FID) 
of the spinP^ In the FID experiment the spin is oriented 
in the xy plane at t = with the field B || z giving the 
major part of the electron precession frequency, and the 



subsequent precession and decay of the average in-plane 
spin is measured. It is crucial to note that the electron 
precession frequency consists of not only the Zeeman fre- 
quency f2, but also the Overhauser field AjJ? arising 
from all the nuclear spins. The current FI D exp eriments 
are either spatial^H or temporal ensemble^^ measure- 
ments. In the first case one addresses many electron spins 
at the same time, each one of them interacting with a 
different nuclear bath. In the second case one measures 
repeatedly the same spin, but the total data acquisition 
time is long enough for the nuclear bath to apprecia- 
bly change its state and probe its phase space. These 
two situations involving spatial and temporal averaging 
are presumably equivalent from a statistical perspective, 
leading to the same characteristic FID time. In both 
cases, the experimental procedure translates into averag- 
ing over the relevant density matrix of the nuclei. The 
rms of the Overhauser field is ~ A/V~N in a thermal 
bath, leading to a strong inhomogeneous broadening of 
the FID signal. 50 The resulting observed decay is Gaus- 
sian, i.e. the in-plane spin decays as exp[— (i/T 2 *) 2 ], with 
a characteristic time-constant of T| ~ y/~N /A being the 
order of 10 ns in III-V quantum dots.£M 1 i3 

The T| decay of FID is a purely classical averaging 
phenomenon, which does not even require the bath to 
have any dynamics. The presence of strong inhomoge- 
neous broadening underlines the necessity of control of 
the random Overhauser shifts of many qubits in a spin- 
based quantum computer architecture, but the actual co- 
herence of a qubit is characterized by a "single-spin" FID 
decay time T^. This is the decay time which one would 
measure for a single spin if the duration of the experiment 
was shorter than the time on which the Overhauser field 
fluctuates appreciably, or for an ensemble of spins hav- 
ing the same Overhauser shift. In the latter case one 
would be dealing with a so-called narrowed state of the 
nuclei. There are a few theoretical prop osals on how to 
drive the nuclei into such a narrowe d state^ 52 * 53 ! and ex- 
perimental progress has been made. 46 48 5 4 ^ In t heoretical 
calculations of FID decoherenc e! 22 1 25 1 26 1 31 1 32 1 33 1 the max- 
imally narrowed states are used (i.e. eigenstates of the 
Overhauser operator), and it is argued that in large bath 
the decoherence dynamics does not depend on the choice 
of the state, with the tacit assumption that the used state 
is a "typical" one (i.e. not one of the few highly polarized 
low-entropy states). In this paper we take an analogous 
approach, and below by FID decoherence we refer to the 
"single-spin" process. 

Whereas the ideas for state narrowing of nuclei are 
very recent, for the past 50 years the preferred way of 
removing the inhomogeneous broadening has been the 
spin echo (SE) experiment in which the electron spin 
is rotated by angle 7r around one of the in-plane axes 
at the midpoint of its evolution. Such a protocol can 
be denoted as t/2 — it — t/2: two free evolution peri- 
ods with a fast external pulse in the middle and readout 
at final time t. It is easy to see that such a procedure 
will remove inhomogeneous broadening (static spread of 
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the precession frequencies). In NMR literature the SE 
decay time is often identified with single-spin T 2 , but 
this holds only when the decoherence due to the actual 
dynamics of the spin and the bath can be treated in 
a Markovian approximation (leading to exp(— tjT^) de- 
cay). This is generally not the case for the spin inter- 
acting with a nuclear bath due to the slowness of the 
bath dynamics and back-action of the electron spin on 
the bathj22HH anc j the Mar kovian decay of spin coherence 
has been predicted 26 ! 33 ! 35 ! only in very specfic conditions 
(long times and large magnetic fields). Thus the char- 
acteristic decay time in the SE experiment, denoted by 
Tse, is most often different from FID T 2 time. The SE 
is the simplest available method of uncovering the non- 
trivial quantum dynamics from underneath the inhomo- 
geneous broadening, and it has been succesfully applied 
in Refs. 1101121 to lateral GaAs quantum dot qubits and 
in Refs. 1 441451 to Si:P qubits (and very recently signif- 
icant progress has been made towards repeating these 
experiments in InGaAs self-assembled dotsj^ where op- 
tical control of electron spin is employed). Consequently, 
a reliable method of calculating Tgg is crucial for estab- 
lishing contact between the experiments and theory. We 
note that the state narrowing protocols of Ref. US] are un- 
likely to lead to a decay time shorter than that obtained 
through the spin echo technique. 

For many years the theories of SE decay due to the 
spectral diffusion (i.e. dipolar interactions) were based 
on stochastic assumptions (see Ref. 122 and references 
therein) Only recently have fully quantum mechanical 
and efficient approaches to the nucle ar bath dynam ics 
under the SE sequence been developedPI 26 ' 27 ' 28 ' 29 ' 30 ' All 
these theories involve some kind of exponential rcsum- 
mation of the perturbation series for the decoherence 
evolution function. The theoretical methods include ex- 
pansion in real-space clusters IHMEU a r elated p air cor- 
relation approximation (PCA) approach] 25 ! 26 * 55 ! and the 
diagrammatic linked-cluster technique! 3 ^ For an excel- 
lent discussion of the relation between these methods 
see Ref. [3H where also a new "cluster-correlation expan- 
sion" is proposed for small nuclear baths, l ike th e one 
interacting with the NV center in diamond. 14 56 These 
theories have shown a very good agreement! 28 ! 30 ! 43 ! with 
the spin echo experiments performed on Si:P system. ^ 4 ! 45 ! 
However, these works have concentrated on a somewhat 
vaguely defined high-field regime, in which the SE decay 
is dominated by spectral diffusion, since the decoherence 
due to the hf- mediated flip-flop of pairs of nuclei is nearly 
completely removed by the SE sequence at B w 10 T. 
This was first observed in numerical simulations,^ and 
then explained in a physically transparent way within the 
PCA approach! 25 ! 26 ! However, this cannot explain the SE 
experiments^ in GaAs at B<0.1 T, underlining the ne- 
cessity for a reliable low-field theoretical approach. 

The calculation of the decoherence in the SE experi- 
ment is in principle qualitatively different than for FID, 
since one is considering the driven evolution of the sys- 
tem. Once the n pulse driving is incorporated into the 



theoretical framework, one can generalize the calculation 
to the evolution driven by sequences of 7r pulses more 
complicated than the SE. The application of such dynam- 
ical decoupling (DD) sequences is expected to prolong the 
coherence timePZEIMl The above-mentioned analytical 
theories were used to calculate this effect for the electron 
spin driven by the classical Carr-Purcell-Meiboom-Gill 
sequence®!, by the concatenated DD ) 61 ! 62 ! and by Uhrig's 
DD sequenced Numerical investigations of the very low- 
field regime (inacce ssible to all the mentioned theories) 
have also been done.HiBH 

Compared to the well-understood status of the spec- 
tral diffusion theory of decoherence due to dipolar in- 
teractions in the bath (relevant at large B), the relation 
between different analytical theories of dephasing due to 
hf interactions alone is less-well established. In the first 
studies of the electron spin dynamics due to hf inter- 
action, the pro cesses o f electron- nuclear spin flip were 
treated directlyj- 22 ! 23 ! 24 ! i.e. the perturbation expansion 
was done in orders of interaction V s { ~ (S + J~ + S~ J + ), 
and in the second order of Generalized Master Equation 
(GME) approach no genuine dephasing (i.e. decay to zero 
of the transverse spin) has been obtained.^ The FID de- 
cay to zero was obtained from V^f using the eq uation s of 
motion method for spin correlation functions, ! 31 ! 32 ! and 
then by the GME approach in the fourth order with 
respect to V^f P^l i n the meantime, a physically intu- 
itive PC A solu tion for both FID and SE decay has been 
proposed] 25 ! 26 ! and later deriv ed using the diagrammatic 
linked cluster expansion ! 30 ! 34 ! In these papers the t^f in- 
teraction was not treated directly, but an effective inter- 
nuclei interaction due to virtual electron-nuclear spin 
flips was used. This hf- mediated interaction approach 
was also used succesfully within the GME frameworkPS 
In all of the above papers, the condition of applicability 
of the theory was either explicitly or implicitly given as 
A/il < 1, which corresponds to B larger than a few Tesla 
in GaAs. At the same time, it has been well established^ 
that the real electron spin-flip is negligible as long as 
5<Cl, which is a much less restrictive criterion. In fact, 
this criterion was suggested to be the relevant one in 
Ref. [25] where PCA method was introduced. Finally it 
should be mentioned that the limit of very low B fields 
corresponding to 8 > 1 (B less than a few mT i n GaAs) is 
accessible only to exact numerical simluations! 39 * 40 ! 41 * 42 ! 

In this paper we present a detailed derivation of a 
theorjEH of pure dephasing under any sequence of ideal it 
pulses, which is based on the effective nuclear Hamilto- 
nian approach in which the electron-nuclear spin flips are 
eliminate d in favo r of hf-mediated interactions between 
the nuclei.^ 33 * 38 ! We argue that for nuclear bath being 
unpolarized and uncorrelated (its density matrix being a 
thermal one) this theory is applicable at experimentally 
relevant low magnetic fields (when A/il > 1), with the 
controlling small parameter S defined in Eq. |2]). This 
result is firmly established only up to the time scale of 
N/A ~ 10 fis (for N= 10 6 ), but at low B most of the de- 
cay occurs on this time scale. The key insight leading to 
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the solution of the decoherence problem is the long-range 
character of the hf-mediated interactions, which couple 
all the N ~ 10 5 — 10 6 nuclear spins appreciably inter- 
acting with the electron spin. This allows us to identify 
the so-called ring diagrams as the leading terms in the 
perturbation expansion of the decoherence function, with 
corrections to them suppressed as l/N. This is analogous 
to the summation of the ring diagrams in the cluster ex- 
pansion of the partition function of the high-density Ising 
model well known from classical statistical mechanics.^ 
A similar approach (termed "the pair correlator approx- 
imation") has also been used recently to calculate the 
influence of the hf-mediated interaction on the electric 
dipole spin resonance in quantum dots.^ Taking into ac- 
count the ring diagrams only, we can sum all the terms 
in the cumulant (linked cluster) expansion of the deco- 
herence function, obtaining closed formulas for decoher- 
ence in FID, SE, and other experiments involving more 
7r pulses. We compar e our FID results with the experi- 
ments on InGaAs dots in which the nuclear state was 
effectively narrowed by driving an ensemble of dots with 
a sequence of optical pulses. We also identify the dom- 
inant processes which lead to SE decoherence at these 
low fields, which are the hf-mediated flip-flops of nuclei 
having different Zeeman energies, i.e. nuclei of different 
isotopes of Ga and As. This crucial role of hetero-nuclear 
character of the bath in III-V materials at low B fields has 
not been appreciated before. At high magnetic fields this 
leads to a small oscillation superimposed on the SE signal 
due to the spectral diffusion, while at low fields it leads 
to a practically complete decay of the SE in less than a 
microsecond in large GaAs dots. This result might pro- 
vide an explanation for the recent measurements of the 
single-spin SE decay on time-scale of a microsecond in 
gated GaAs quantum dotsP^l 

Our work, developed in details in this article, brings 
the theory of quantum dephasing of the localized electron 
spin qubit due to its hf interaction with the nuclear spin 
to the same level of sophistication and quantitative depth 
that the current theories of spectral diffusion due to the 
dipolar nuclear flip-flops have. As such, we can carry out 
quantitative comparison with experiments performed at 
low magnetic fields, where the hf coupling would domi- 
nate over the spectral diffusion process. It is therefore 
gratifying that our theory gives reasonable quantitative 
agreement with the existing experimental data on spin 
dephasing down to a field as low as 10 mT. More exper- 
imental data would be required for any further improve- 
ment in the theory. 

The paper is organized in the following way. In Sec. |TT] 
we describe the full Hamiltonian of the electron spin and 
the nuclear bath, and derive the effective pure dephas- 
ing Hamiltonian containing a hierarchy of multi-spin hf- 
mediated interactions. In Sec. |III| we present a general 
formal framework for calculation of pure dephasing un- 
der the action of a sequence of ideal tt pulses. The main 
theoretical development is presented in Sec. |IV| where 
we identify the most important terms (corresponding to 
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TABLE I: Typical values of energies in the Hamiltonian in 
units of angular frequency (ft=l) and temperature (fcg = 1). 
The Zeeman energies Q, and u>i are given for B — 1 T. The 
dipolar interaction Bij is given for two nuclei separated by 1 
nm distance. GaAs parameters are taken from Ref. IBB] (see 
also Ref. I20|l . The values of A a for In are taken from Ref. 1261 



the ring diagrams) in the perturbation series for the de- 
coherence function in the case of long-range interaction 
between the nuclear spins, and we show how they can 
be resummed. This approach is then used in Sees. [V] 
and |VI| to calculate the decoherence in single-spin (nar- 
rowed state) FID and SE experiments, respectively. Our 
FID calculation is compared with other recent theories 
in Sec. |VB| and in Sec. |V C| we show that the ring di- 
agram approach to FID agrees with a result obtained 
in a completely different fashion in the special case of 
uniform coupling of the nuclei to the central spin. In 
Sec. |V D] we present our results for FID in GaAs and In- 
GaAs quantum dots. Sec. |VI| contains the theory and 
results for SE decay in low B fields, when the interaction 
between the nuclei of different species gives the dominant 
contribution to SE decoherence. The theoretical predic- 
tions for multi-pulse dynamical decoupling sequences in 

In 



the same magnetic field regime are given in Sec. VII 
Sec. IVIIII we discuss the corrections to the lowest-order 
hf-mediated interaction and argue that our approach is 
well controlled on experimentally relevant time-scale as 
long as S < 1 (B » 3 mT in GaAs dots with N ~ 10 6 
nuclei) . 

Before proceeding further, let us establish the nomen- 
clature. By the short-time regime will mean t <C N/A 
(which corresponds to t <C 20 /is in GaAs with TV = 10 6 
nuclei), while t^$>N/A is the long-time regime. By high 
magnetic field B we will mean B such that A/ ft < 1 
(corresponding to B > 3 T in large GaAs dots), while 
low B will mean that A/Cl > 1, but 6 = A/Cly/N < 1. 
The regime of very low B corresponding to 5 > 1 remains 
outside of the realm of validity of the theory presented 
below. 
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II. THE MODEL 

The full Hamiltonian of the electron spin interacting 
with the the nuclear bath is given by 



H — Hz + -ffdip + H] 



hi 



(3) 



with the terms corresponding to Zeeman, dipolar, and 
hf interactions, respectively. In order to describe the nu- 
clear system in III-V semiconductors we have to account 
for the existence of nuclei of different species (isotopes) , 
which we label by a. To the lattice sites i we assign 
randomly the nuclear species a[i], and in the subsequent 
calculations the average over possible nuclei on a given 
site will be taken implicitly, using the weights n a denot- 
ing the number of a species nuclei per unit cell. Nuclei of 
different a have different gyromagnetic factors, resulting 
in different Zeeman splittings, hf couplings, and dipolar 
couplings. Thus we write Zeeman and dipolar energies 
as 

H z = QS Z + J2^[i\Ji , (4) 

i 

ffdip = J^hjiJfjr -2JfJ/) , (5) 

and the hf interaction is given in Eq. (JT|) . In these equa- 
tions fl is the electron's spin Zeeman splitting, u) a u] are 
nuclear Zeeman splittings for nuclei at the i-th site, and 
bij is the secular dipolar interaction between the nuclei, 
with a[i] — a[j] implied.^ The secular approximation 
to the dipolar Hamiltonian is justified by the fact that 
bij <C u)i,u)j and bij <C uj a p = to a — ujp for a ^ f3 for 
magnetic fields of interest here (see Table [I]) . The hf cou- 
plings Ai in Eq. (JT|) are proportional to the total hf in- 
teraction energy A a [i] of the nuclear species a[i] at site 
i, and the square of the envelope function of the electron 



(6) 



with the envelope normalization given by: 

f \nr)\ 2 dr = v O $> = n C) (7) 
Jv 

where V is the volume of the quantum dot, vq is the 
volume of the primitive unit cell (^o = Oq/4 with ciq be- 
ing the lattice constant in zinc-blende or diamond struc- 
trure), and n c is the number of nuclei in the primitive 
unit cell (equal to 2 in GaAs) . In terms of microscopic 
parameters, we havtP 



Ai = ^^ 2 lSlJa\u a \ 2 



(8) 



with [Xq the vacuum permeability, 7s and jj a are the 
electron and nuclear spin gyromagnetic factors, respec- 
tively, and u a is the amplitude of the periodic part of the 
Bloch function at the position of the nucleus of a species 



(the normalization is J |w(r)| 2 dr = l). The values of A a 
for Ga, As, and In are given in Table [i] 

It will be convenient later to use the following defini- 
tion of N, the number of nuclei interacting appreciably 
with the electron: 



N 



Ei/i = /l*(r)| 2 rfr 
/|*(r)|4dr 



(9) 



from which we get that when summing over the nuclei of 
one species a we have 



A 



*A 2 JN 



(10) 



where we have used the fact that J2ie a can be replaced 
by (n a /n c ) when dealing with a large dot (i.e. when 
there are enough nuclei for the distribution of Ai to be 
self-averaging). In the same way we get the total hf in- 
teraction energy 



A — ^ Ai ^ fi a A a , 



(11) 



which, using the parameters from Table [T] gives A = 
83 /ieV in GaAs. Using the effective g-factor of g e s = 
—0.44 in GaAs we get that A/SI = 1 for magnetic field 
B = 3.25 T. On the other hand, the rms of the Over- 
hauser field in an unpolarized thermal ensemble of nuclei 
is A/vN, which for N — 10 6 corresponds to magnetic 
field of about 3 mT. Thus, in the following by high fields 
we will mean B > 3 T and by low fields we mean B < 3 
T (but larger than ~ 10 mT so that S is indeed a small 
parameter) . 



A. Distribution of hyperfine couplings Ai 

In the following we will assume that the envelope wave- 
function of the electron in the dot is 



*(r) 



2^ 



ttL 2 



exp 



y 



2L 2 



(12) 



with the thickness of the dot zq and the Fock-Darwin 
radius L. For this choice of the envelope, from Eq. (|J 
we get that N = Azq-kL 2 /3i/o- This number is of the same 
order of magnitude as N m — zqttL 2 /2i/q defined by 



max Ai = A/N n 



(13) 



so that whenever we encounter a sharp inequality involv- 
ing Ai, we can replace Ai in it with A/N . 

For systems with a large number of nuclei we can make 
a continuum approximation for the distribution of A^. 
We write for any function F(A{) 



V F(AA = n a [°° F(A) Pa (A)dA 



(14) 
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where the "density of states" of the hf couplings is 

p a (A) = — [ S(A-A a \*(r)\ 2 )d 3 r , (15) 
^0 Jv 

so that p(A)dA is the number of unit cells in which the 
hf interaction energy is in the range [A, A + dA\. For the 
electron wave-function from Eq. (1 1 21) we get 



development of a theory of decoh erence d ue to various 



3 N 3NA r^f&A a . , 

p Q (A) = --arccos^— 6^--A) , (16) 

where &(x) is the Heaviside step function. d 



Effective Hamiitonian 



We write the hf interaction as 

H u = J2AiJfS z + V si 



(17) 



where = gX^A^S+Jr + S'jf) is the spin-flip 
part. In a finite magnetic field this interaction was 
showrP^ 38 ! 39 ! to lead to relaxation of the S z component 
of the electron spin by a small quantity S 2 = A 2 /Nil 2 . 
In a large magnetic field when £l^> A, the direct contri- 
bution from t^f can be neglected when considering the 
electron spin dephasingP^ However, while the real S — J 
spin flip is forbidden, higher order virtual processes in 
which the electron spin is flipped multiple times while 
having the same initial and final state, are still allowed. 
The existence of such processes is clearly brought to 
light by performing a canonical transformation remov- 
ing V s f from the Hamiitonian. Such a Schrieffer- Wolff 
type transformation leads to an appearance of effective 
hf-mediated interactions between the nucleiPSE 33 | 35 | 3S | 69 1 
We give a more detailed derivation of this transformation 
in Appendix [A] Here let us mention that in deriving the 
effective interactions we have assumed that f2 is much 
larger than the rms Overhauser shift in the nuclear state 
under considration. Since we are interested in thermal 
and unpolarized states of nuclei, this means f2^> A/vN 
(i.e. <5<§;1). Under this transformation the Hamiitonian 

changes to H = e~ s He s and the states are transformed 
to \<f>) = e~ s \cf>) with a carefully chosen unitary opera- 
tor e~ s . An important observation is the fact that the 
transformation of the states can be showrP^l to lead to 
a partial decay (visibility loss) of the coherence occur- 
ing on time-scale t < N / A and having the amplitude of 
6 2 , in qualitative agreement with results from the GME 
theoryP2 This leads to the conjecture that the "effective 
Hamiitonian" approach is physically sensible at least up 
to some time-scale as long as 5 <C 1. We elaborate on 
this issue in Sec. | VIII] where we show that at the experi- 
mentally relevanlP^ 1 time-scale (less than ~ 10 /j,s in GaAs 
with N— 10 6 nuclei) this condition is in fact sufficient for 
our approach to be well-controlled. Since this discussion 
involves time-dependence, it has to be preceded by the 



mfvr 



terms in H, which is given in Sees. 

In all of the following we will assume 8 <C 1 and neglect 
the transformation of states accompanying the transfor- 
mation of H. Using the method described in Appendix 
[A] we arrive at the effective pure dephasing Hamiitonian 



H = H Z + J2 A iS*Jf + Hdi P + E ^ 



(n) 



(18) 



n=2 



where represent nuclear interaction terms generated 
in n-th order with respect to V 8 f. The terms can 

contain only S z operator of the electron spin, and thus 
in Eq. ( 18 1 we have an effective pure dephasing Hamilto- 



In the seco nd order we get a previously known 
resulPIMm 



#(2) 



EA 2 „ 
4ft Jl 



(19) 



The first term is a very small correction to Zeeman en- 
ergy of nuclei, corresponding to the renormalization oji 
in the original Hamiitonian by Aw| 2 ' = — Af/4£l, which 
is a negligible correction for B fields above a millitesla. 
The second term is also rather small, and for J — 1/2 it 
amounts to a renormalizatiorP^ of electron spin splitting: 
AO (2) w A 2 /4:Ntt. The third term is the hf-mediated 
two-spin (2s) flip-flop interaction, which will be the fo- 
cus of most of our attention: 



H. 



(2) 



r r 



(20) 



where Bij = AiAj/AQ. Since the nuclear interaction in 
this term is multiplied by electron S z operator, we will 
call terms of this kind the ^-conditioned ones. The cru- 
cial characteristic of is that it describes a non-local 
interaction: we have of the order of iV 2 pairs of nuclei 
coupled by appreciable interaction constants ~ A 2 /N 2 il. 
Thus, even if the typical interactions energy of any given 
pair of nuclei is small, one has to remember that this 
weakness of couplings can be overcome by the number of 
coupled pairs of nuclei. 

The full expression for purely hvperfinc contribution 
to H^ 3 ' is reproduced in Appendix \A\ In this order we 
also obtain terms due to mixing of hf with dipolar inter- 
actions and nuclear Zeeman terms. These turn out to be 
small corrections to the pre-existing terms or interactions 
derived from V s { only. Involving the dipolar interactions 
makes also the effective coupling more local, while the 
most important interactions are the "maximally nonlo- 
cal" ones, such as the three-spin (3s) term 



E 



AiAjAk 

2n 2 



J A J A J h 



(21) 
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Another important feature in the third order is the ap- 
pearance of the S^-independent hf interaction: 



H, 



(3) 



A, 



im 2 



(22) 



(2) 



The 5 z -conditioned two-spin interactions in arc 
small corrections to the flip-flop interaction from H. 
since their coupling constants are scaled down by 

In the following Sections we will concentrate on the 

~ (2) 

S^-conditioned two- spin interaction from H^J considered 
previously in Rcfs. 25 26 33 61, and on ^-independent 

~ (3) 

H-2s terms. The influence of the higher-order multi-spin 



terms (such as H^ s ') will be discussed in Sec. VIII 



III. GENERAL APPROACH TO 
DECOHERENCE UNDER PULSE SEQUENCES 

We write the pure dephasing effective Hamiltonian as 

H = H + V = H Zn + 2S z H A + 2S z V 1 + V 2 , (23) 

in which H Zn is the Zeeman energy of nuclei, 2S Z Ha is 
their diagonal hf coupling to the electron spin (Ha = 
J2i ^t^Di and all the interactions between the nuclei 
are divided into S^-conditioned ones in V\ and S z - 
independent in Vz (which could include also the dipo- 
lar interactions). The electron Zeeman energy has been 
dropped from H since it contributes only a simple phase 
factor of Sit to the FID signal and completely cancels 
out from the expressions for SE decay. We consider the 
evolution operator U(t) of the whole system under the 
action of a series of ideal (<5-shaped) tt pulses applied to 
the electron spin. For pulses corresponding to rotations 
by angle tt about the x axis we have 



U(t) = {-i) n e 



-iHi 



-iHi 



(24) 



with n being the number of applied pulses, ti being time 
delays between the pulses, and the total evolution time 
t = y"ffj Tf. For example, the spin echo corresponds to 
?i=l and T\ = T2 = t/2. 

The coherence of the electron spin is given by the off- 
diagonal element of its reduced density matrix, given by 



p + -(t)=(+\Tr B p(t)\-) 



(25) 



where |±) are the eigenstates of S z , p is the density 
matrix of the whole system (clcctron+nuclei), and Tre 
denotes tracing over the bath degrees of freedom. We 
assume the factorizable initial density matrix given by 
/5(0) = Ps( 0)pb(0) - With this we arrive at the following 
expressio n 1 60 * 61 ! 70 * for the decoherence function W(t) for 
even n: 



W(t) 



P+-(t) 
P+-(0) 



Ul(t)U+(t) 



(26) 



and for odd n there is p (_(0) in the denominator. This 

definition ensures W(0) = 1. In Eq. (26 1 we have (■••) = 



Ttb{pb(0)..-} denoting the average over the states of the 
nuclei, and U± are evolution operators of the nuclear sys- 
tem conditioned on the state of the electron spinP^l 

U(t)\±) = \±)U±(t) for n even, (27) 

U(t)\±) = |=F>&f(*) for n odd. (28) 

For example, we have for FID 



e iff -V 



iH+t 



(29) 



with H±—Hzn ± Ha ± V\ + V2, and for SE we have 

fjSE _ e -iH±t/2 e -iH T t/2 



(30) 



and the expression for Wse(^) is obtained by plugging 
these into Eq. (26 1. 

Note that in Rcf. 35] we defined W(t)=2\p + -(t)\. As 
we discuss below, the two definitions are equivalent for SE 
and other dynamical decoupling sequences, but for FID 
the earlier choice meant discarding of the phase dynam- 
ics of /9_| , which will turn out to be nontrivial. Experi- 
mentally, the current definition corresponds to measuring 
(S x {t)) or (S v (t)), while the previous one corresponded 
to measurement of (S x (t)) 2 + (S v (t)) 2 . The latter would 
require averaging over a series of measurements in which 
cither x or y components of the spin are measured. 



A. Decoherence function on a contour 

The evolution under the influence of pulses can be re- 
expressed by making Ha and V\ time-dependent. Let us 
define the filter function in the time domain f(t;r): 

n 

f(t; t) ee (-1)" ^(-l) fe 6(r - t k )Q(t k+1 - t) , (31) 

k=0 

where t k with k—l...n are the times at which the pulses 
are applied, io = and t n+1 =t, so that f(t; r) is nonzero 
only for r € [0,t1- Any pulse sequence can be encoded by 
this function j 71 * 72 le.g. the CPMG sequence^ corresponds 
to t/. = t(k — 1/2) jn. Using this function we can write 
the closed formula for U± as 



U±{t) = T exp 



1 I (H Zn + V 2 ±f(t;T)[H A 
'0 



V^dr 



(32) 

where T is the time ordering operator. Then we trans- 
form into the interaction picture with respect to V(t) = 
V% ± fit; tf)Vi. The evolution operator in the interaction 
picture, S±(t), is defined by 



U±(t) = exp 

EE U£(t)S±{t) 



i I Hq (t) d,T 




S±(t) , (33) 
(34) 
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where H^(t!) — Hz n ± f(t;r)HA, so that the time- 
ordering is actually unnecessary and in expression for 
U^{t) and it is given by 



U±{t) = exp 



-iH Zn tTiH A / /(i;r)dr 
Jo 



From this we get 

S±(t) = T exp 
Sl(t) = Texp 



-i / V ± (r)dT 
Jo 

i [ V ± {r)dT 
Jo 



(35) 

(36) 
(37) 

(38) 



Then, for any balanced pulse sequence (defined by 
J f(t;r)dT — 0, which does not hold for FID, see 
Sec. IIICl, the decoherence function becomes: 



where T is the time anti-ordering operator and 
V±(r)=(u ± (r)yVU ± (r). 



W(t) = (sl(t)s + (t) 



(39) 



In Eqs. (26 1 and (39) one can see that W(t) is the 



average of the evolution of the bath first from to time 
t, and then back (under a different Hamiltonian) from t 
to 0. This structure calls for use of a closed-time-loop 
technique well-established for nonequilibrium problems 
in many-body theory.^ Such a Keldysh-like str ucture has 
been used previously in decoherence problems! 30 * 34 * 70 ^ 
Introducing the notion of ordering the operators on the 
Keldysh contour (see Fig. [IJ allows us to write W(t) in 
a more compact way: 



W(t) = (T c exp 



i I V{T c )di 



(40) 



In the above equation Tc denotes the contour-ordering 
of operators,^ and r c = (r, c) is the time variable on the 
contour, with c = +(— ) on the upper (lower) branch of the 
contour and r being the corresponding real time value. 
The contour integration of a function F(t c ) — F(t, c) is 
given by 

p />oo />0 

/ F(r,c)dr c = / F(r,l)dr+ / F(r,-l)dr . (41) 

JC JO Joe 

In terms of the free-evolution operator on the contour: 



^o(tc) = exp 



iH Zn r - ic / f(t; t')H A dt' 



(42) 



we have the interaction operators 

Vi(r c ) = c/(i;r)$(T c )ViWo(T c ) , (43) 
V 2 (r c ) = uI{t c )V 2 Uo(.Tc) ■ (44) 



c= + 



T = 



r = t 



b) 



-1... 



c) 



f(t;r) for SE 



f(t\ r) for 2-pulse CPMG 



FIG. 1: a) The closed loop contour along which the operators 
V(t, c) are ordered, b) The plot of the time-domain filter 
function f(t; r) for the Spin Echo sequence, c) The same for 
2-pulse CPMG sequence. 



Note that the sign of the Vi interaction is opposite on two 
branches of the contour. The Vi 2 interactions are given 
by the same formulas as the Schroedringer picture Vi,2> 
only with the spin operators replaced by their branch- 
and time-dependent counterparts: 



Jf{r c ) = Jf exp 
and J z (t c ) = J z . 



A ■ 

±iLU 3 T±ic I f(t;t')^-dt' 
2 



,(45) 



B. Averaging over the nuclear ensemble 

The averaging over the nuclear ensemble is performed 
here using a thermal density matrix of nuclei pb — 
e~P HB /Z, where Z is the partition function. Due to 
the smallncss of the energy scale associated with nuclear 
Hamiltonian Hq = Hzn + -Hdip, at the experimentally 
achievable temperatures we approximate it by pB ~ 1/Z. 
Thus, we treat the case of the "infinite temperature" un- 
polarized nuclear bath. Then, the average in Eq. ( |39[ ) for 
balanced pulse sequences is a real number and there are 
only even orders of V in the perturbation expansion of 
W(t) (see Ref.HD. 

An important feature of the thermal nuclear bath is 
that it is uncorrelated, i.e. the average of a product of 
nuclear operators corresponding to disjoint sets of nuclei 
factorizes: 



) J/ 

J, 



,J 



,)> = 



(46) 



where i ni } and {ji,--,jn 2 } are two non- 

overlapping sets of nuclei. Below we will extensively use 
this factorization in the case when among the averaged 
operators there are only two ( and J~ ) corresponding 
to a certain nuclear spin i. In such a case we have 



(...J±...J?...) = (J+jr)(...) 



(47) 
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where 



(J^JD = -Ji(J i + l) = a i , 



(48) 



for a nuclear spin J^. In the most relevant case here of 
GaAs quantum dots, we have Ji — 3/2 and en = 5/2 for 
all the nuclei, while for nuclei of Indium we have Ji~ 9/2 
and Qi = 33/2. 



C. Decoherence function for Free Induction Decay 



In the case of FID, the diagonal hf interaction Ha 
in Un does not cancel out from W{t), and instead of 



Eq. (39) we obtain 



W Fm (t) = (SL(t)e 



-2iHAt 



s+(t) 



(49) 



Using the "infinite temperature" nuclear thermal ensem- 
ble we arrive atP^ 

j 

where \J) = ®, \jf ) are the product states of the eigen- 
vectors of Jf operators of all the nuclear spins, and £1 j 
is the Overhauser shift of the electron frequency for the 
nuclear state 



Slj = {J\2H A \J) = J2 A i3 



(51) 



The form of the decoherence function in Eq. ( 50 ) is quite 
different from the one in Eq. (39 1. However, it has been 
shown in Refs. [231 and I2U1 that the matrix element in this 
equation for practical purposes does not depend on the 
state \J), at least as long as we consider the "typical" 
states with polarization close to zero. This follows from 
the fact that there are many nuclei in the bath making 
it self-averaging. Let us note that this is not the case for 
dilute baths, such as the dipolarly coupled spin bath in- 
teracting with the NV center in diamond, where different 
spatial realization of the bath can lead to very different 
FID signals P 

Consequently, in the dense hf-coupled bath we make an 
approximation of replacing the (J7"| ... \J) matrix element 
with its ensemble average, obtaining 



w FlD (t) » j2> 



—iQ ft 



x (SL(t)S+(t) 
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Wfid(*) x W Fm 



(52) 
(53) 



In the above Equation, W FIU gives t he fast FID de- 
cay due to inhomogeneous broadening! 23 1 2615 ^ W Fm ~ 
exp[-(i/T|) 2 ] where T 2 * w y/NA~ 10 ns, while W| ID is 
the decoherence function for a single spin. The latter 
quantity, which has the same form as Eq. ( 39 ) , is of main 
interest here. 



The measurement of the single spin FID is in princi- 
ple possible by devising an experimental procedure which 
leads to "nar rowing" of the distribution of the nuclear 
states IHl 52 l 53 l with only states having the same Clj con- 
tributing in the fully narrowed situation. There has re- 
cently been experim ental progress in performing such 
experiments l 46 l 48 EH with an estimate of the single-spin 
Ti time of about 3 fjs in InGaAs dots.™ 

The above-described approach for calculation of W F1D 
is equivalent to calculating the FID evolution for a single 
"typical" \J) state, which was employed in Refs. |2"21 and 
33 While in most of the following we will concentrate 
on the calculation in which W F1D is calculated by full 
ensemble average, with Slj shifts factored out according 
to Eq. (53 1, in Sec. VC we will compare this approach 



with another method in which we explicitly worked with 
a narrowed nuclear state. 



IV. RING DIAGRAM EXPANSION OF 
DECOHERENCE FUNCTION 

The theoretical framework presented so far is very gen- 
eral. The calculation of the decoherence function given 



in Eq. ( 39 1 has been achieved for the nuclei interact- 



ing via dipolar interactions using different versions of 
expon ential resumm ation of the perturbation series for 
W(t) )2512 6 | 27 | 28 | 29 | 3 0] s j m ii ar approach has been used 
for the case of nuclei interacting via the low est-order S z - 
conditioned hf-mediated interaction * 25 ! 26 * 3 ^ but only the 
lowest-order terms in the cumulant expansion have been 
retained in these calculations. In this Section we show 
how to resum all the leading (in terms of 1/N expansion) 
terms in the perturbation series for W(t). We start with 
an illustrative example calculation of the lowest-order 
terms in Sec. |IV A| give the full theory of the resum- 
mation of ring diagrams in Sec. |IVB[ discuss its relation 
to previous works in Sec. IV C and finally present a nu- 



merically efficient method of calculating the sum of all 
the ring diagrams in Sec. |IVD| 

In this Section and in Sees. fVl and IVII we will concen- 
trate on two-spin hf-mediated interactions, such as the 
leading ^-conditioned interaction in H given in Eq. ( 20 1 



and the ^-independent term from Eq. (22 1. The role of 



multi-spin interactions will be discussed in Sec. VIII 



A. Lowest order contributions to W(t) 

For concrctcness, let us use the two-spin interaction 
from Eq. ( 20 ) . The corresponding interaction on the con- 



tour is given by 



V(r c ) = cfit-T^B^Jt^J;^) , (54) 
= Y, V MJ+j- , (55) 
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with the contour-time dependence of the ladder opera- 
tors given in Eq. (45) and with B^j = AiAj/4£l. For the 



two-spin interactionfrom Eq. ( 22 I the calculation is anal- 



ogous, only there is no cf(t; t) factor in Eq. (54 1 and the 
couplings Bij are different. 

Perturbation expansion of W(t) defined in Eq. (40) 
leads to the series W = 1 + £ fc=2 W W , with ~ (W). 
In fact, for SE and any balanced sequence of pulses only 
the even orders are non-zero and only for FID we need 
to consider odd k. In order for the product of V not 
to average to zero, all the raising operators J 2 + have to 
be paired up with lowering operators jr. In the second 
order we have 



1 



c Jc 



dr Cl / dT C2 (T c V(T Cl )V(r C2 )) , (56) 



where we used the notation r Cj = (rj,c,*). The averaged 
expression is 

(T c V(t Ci )V(t C2 )) = E V fcI (r Cl )V mn (r Ca ) 

k^l rriy^n 

(T C J+J-J+J-) , (57) 

and since k ^ I and m=/=n there is only one way of con- 
tracting the operators: 



( J k J l J m J n ) - SknSlm{J k Ji J k ) 



(58) 



Let us note now that due to k ^ I constraint, the ladder 
operators effectively commute under the sign of the av- 
erage: the only non-zero commutator is [jf , jr] = 2 jf 
with i = k,l, and while moving the operators around we 
can create at most one j z operator. The average of an 
expression with a single j 2 being the only operator per- 
taining to the i-th spin can be factored as in Eq. ( 46 ) 



and one can see then that the whole expression is equal 
to zero since (jf) = in the unpolarized bath. Therefore 
the spin operators are effectively commuting under the 
average, and we can remove the time-ordering arriving 
at 



(59) 



k,l 



where we have defined the matrix T)y(t): 



T kl (t) 



(1 - S k i)^a k ai / d,T c V kl (T c ) , (60) 
Jc 



with etfc defined in Eq. ( 48 ) and Vki defined in Eq. ( 55 1 



This quantity can be readily evaluated for any two-spin 
interaction and any pulse sequence. 

An expression from Eq. ( 59 1 can be represented graph- 
ically by a diagram shown in Fig. |2^i, where the nuclei 
are represented by the dots, and the interactions inte- 
grated over the time along the contour (i.e. the elements 
of the T k i matrix) are represented by lines. In the cor- 
responding analytical expression each dot contributes a 



sum over all the N nuclei, with the restriction that only 
terms with all the indices corresponding to distinct nuclei 
are summed over in the whole expression. Physically the 
diagram in Fig. [2^i describes a process in which a pair of 
nuclei flips twice during their evolution along the "time" 
contour and comes back to the original state (so that the 
whole process survives the averaging procedure). Simi- 
larly, the higher order diagrams in Fig. [2] correspond to 
processes of propagation of a nuclear spin flip along a 
closed loop involving more than two nuclei. While the 
probability of a single such event quickly decreases with 
the number k of the nuclei involved, the number of such 
events grows as N k . The ring diagrams represent the 
sums over all such events. 

In an analogous manner one can calculate the third 
order term W^ 3 \t): there we have two possible contrac- 
tions, but similarly to the case of we have then a 
sum over three distinct nuclear indices, and the ladder 
operators are again effectively commuting. The final re- 
sult is 



{-if 
3! 



2 T k i(t)T lm (t)T mk (t) , (61) 

k^l^m 



which is represented graphically in Fig. |2]d. 

However, the higher orders of perturbative expansion 
of W(t) do not have the same simple structure. Below we 
discuss the approximation under which all the most im- 
portant contributions to (t) can be easily calculated 
in any order k. 



B. Exponential resummation of W(t) 

When calculated exactly, the higher order terms in ex- 
pansion of W(t) are quit e cumbersome, since the dia- 
grammatic expansion^ 30 ! 34 ! is complicated by the lack of 
simple Wick's theorem for spin operators, and in order to 
account for all the terms one has to introduce more com- 
plicated diagram construction rules than the ones given 
for diagrams corresponding to and M^ 3 '. Now we 

will describe an approximate way of resumming of the 
perturbation series for decoherence due to the non-local 
interactions between the nuclei. 

The term has the following structure: 



(62) 



so that it is a sum of [N(N— l)) k averages of time-ordered 
products of operators. These averages are non-zero only 
when all the J + and J~ operators are paired up, so that 
there are at most ~ N k terms contributing to the sum. 
From these sums over k nuclear indices we now take only 
the terms in which none of the indices is repeated. As- 
suming that all the terms are of the same order of mag- 
nitude, this amounts to at most 1/N error compared to 
the exact expression. The simplification following from 
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this step is tremendous: now, as in the cases considered 



in Sec. IV A the operators effectively commute under the 
average and we can get rid of the time-ordering. Let us 
stress that in this approximation we arc neglecting all the 
terms which correspond to diagrams having more inter- 
action lines than distinct nuclei. Rules for constructing 
these subleadin g (for long-range interaction) diagrams 
are complicated j 3 ° l 34 l but fortunately we do not need to 
use them here. 

Now, only the way in which the spin operators are 
paired up matters. In the 4th order we have two distinct 
patterns of contractions, see Fig. [2]:. There are 6 ways of 
contracting the operators in such a way that the interac- 
tions Bij are connecting all the nuclei in a closed cycle 
(a four-nuclei ring), and 3 ways of forming two disjoint 
cycles (2-nuclci rings) . The sum is over 4 distinct nuclear 
indices, but we will invoke the 1/N expansion again and 
allow for the indices to repeat if they appear in disjoint 
rings connected by interaction, leading us to 



4! 



(6R 4 + 3Rl) 



(63) 



and in a similar way we get for the higher order terms: 

(64) 



W i5} (t) « ^^(4LR 5 + 20^2) 



W {6) (t) 



5! 

H) e 

6! 



(51^6 + 90^4^2 + 5!!^) , etc (65) 



where we have defined the ring diagrams R k : 



Rkit) 



E 



Ti x i 2 (t) ...Ti k i 1 (t) 



TrT(t) k . (66) 



The sum on the left-hand side of the above expression is 
restricted to distinct nuclei (i.e. i n + 1 i m for each n and 
to). Consequently, the approximation in the rightmost 
part of the above formula again introduces a 1 /N error. 

The combinatorics of generating ring-diagram approx- 
imations to higher-order terms in expansion of W(t) 
quickly becomes tedious. For example, in W^ 2k ^ the 
number of contractions leading to R k terms is given by 
(2fc — 1)!!, and the number of Rk terms in is {k— 1)!. 
However, for further progress we only need to know the 
coefficient in front of Rk contribution to W^ k \ since this 
is the only connected term: all the others are products of 
multiple rings. 

Here we use a very g eneral linked-cluster (or cumu- 
lant expansion) theorem ) 77 ! 78 ! according to which for a 
quantity W(t) being an average of a generalized exponent 
(e.g. an ordered exponential of operators) , the logarithm 
of W(t) is given by the sum of all the terms in expansion 
of W(t) which are irreducible, i.e. cannot be written as 
products of averages. In the diagrammatic language this 
means that In W(t) is the sum of all the linked diagrams 
with combinatorial prefactors with which they appear in 
expansion of W(t). A short proof of this theorem (along 
the lines of Ref . 175]) is given in Appendix |b| In k-th order 
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FIG. 2: ( a) a nd (b) Graphical repr esen tation of exact W^ 2 ' 
from Eq. (59 1 and from Eq. (|6l| . (c) Ring diagram 

approximations for from Eq. B63p . (d) Linked cluster 

resummation of W(t). 



of expansion we have then the linked term: 



W ( k ) 

linked jL| 



(-i) k (-i) k 

I (k- l)\R k = ^-R k (67) 



leading to the closed expression for the decoherencc func- 
tion 



W(t) » exp [J2 



-i) k 



R k {t) 



(68) 



\k=2 



where the ~ sign reminds the fact that we have used 
1/N approximation, and that we have replaced the up- 
per summation limit (which is of the order of N) by infin- 
ity. The ring diagram terms R k (t) are given by Eq. (66 1. 



The graphical representation of this equation is shown in 
Fig. [2]i. For SE and other balanced pulse sequences we 
then have 



/ (-l) k 
W cvcn (t) « exp £ y —^R 2k (t) 

\k=l / 



(69) 



These equations arc the central formal result of this pa- 
per. 



C. Relation to previous cluster-expansion theories 
of electron spin dephasing 

Before moving on to discussing efficient ways of eval- 
uating Eq. ( 68 1 , let us outline the relation of this the- 



ory to the approaches previously used to calculate the 
pure dephasing due to the nuclear bath. Retaining only 
the i?2 term in the exponent corresponds to the PCA 
of Refs. |25l26] Note that since W^ = -(1/2)R 2 exactly, 
without invoking 1/N expansion, R 2 is also the first term 
in the linked-cluster expansion for the dipolar flip-flop in- 
teraction between the nuclei, and in Refs. 25 26 61 it was 
used to calculate W(t) for both hf-mediated and dipolar 
interactions. It was noted in Ref. [STJlthat for the long- 
range interaction the diagrams connecting cyclically the 
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maximal number of the nuclei in a given order of expan- 
sion are the most important. These are the ring diagrams 
considered here. 



D. 



Calculation of all the ring diagrams using 
T-matrices 



Furthermore, we can relate Eq. (68) to the real space 



Eq. ( 66 ) can be rewritten using the eigenvalues A; (t) of 



the T-matrix as R k {t) = £f Af (t). Then, from Eq. ([69 1 
we get 



cluster expansion of Refs. [27 28 29 by noticing that R 2k 
are related to the cluster contributions 
Ref.[2E] 



2k 



defined in 



N 



^Wcvc„ = n cx p hE 



If (-1) 



k + 1 



2 ^ k 

k=l 



21, 



(71) 



(-1)* 
2k 



„ N , s N 

R2k = }^v' 2k in the 2k-th order in V, (70) = JJ exp (^_^ m (j + A; (<) 2 ) j = [] JL _ . , ( 72 ) 



where the sum is over all the clusters having 2fc nuclei. 
From the perspective of real-space cluster expansion of 



Refs. 27 28 29 the exponentiation in Eq. ( 68 ) serves to ac 



count, with correct combinatorial factors, for all products 
of contributions from disjoint clusters that must be in- 
cluded in W(t); the 1/nl factor of the exponential expan- 
sion serves to compensate for the n! permutations over 
the n cluster contributions in each product. Errors in this 
approximation arise from extraneous products among 
overlapping (non-disjoint) clusters. Here, these corre- 
spond to the aforementioned 0(1/N) errors. In these 
previous papers the cluster contributions v' 2k were calcu- 
lated for nuclei coupled by dipolar interactions, which are 
not long-ranged in the meaning used here, i.e. they do not 
couple equally all the N nuclei. Thus R 2k with k> 1 are 
not expected to be as good approximations to the exact 
cumulant of order 2k as in the case of the hf- mediated 
interactions. Indeed, for dipolar interactions the above- 
mentioned cluster overlap corrections are important for 
accuracy of higher-order terms in the exponent!^- These 
corrections are related to more complicated diagrams ap- 
pearing in the linked-cluster expansion and recently an 
approach in which they are avoid ed ha s been proposed 
for the case of small nuclear bathsJ^EH 

In the case of dipolar interactions the contribution 
of R 2 w &s gi vm g a. certain decoherence time-scale, and 
the contributions of larger clusters were negligible at this 
time-scale. This can be traced back to the fact that the 
number of relevant terms, which involve close-by nuclei 
with large couplings, was of the same order of N in all 
considered orders of expansion. In such a case the cor- 
rect small parameter for the expansion is bt, with b being 
the largest possible dipolar coupling, and consequently it 
is enough to consider only small clusters of a few nuclei. 
The situation is different for hf-mediated interaction: the 
number of terms to consider in fc-th order is N k , and at 
low magnetic fields (when the hf-mediated couplings are 
not very small) it is not a priori clear that one can con- 
sider only the smallest clusters of nuclei. In fact, in the 
following we will show that at experimentally relevant 
magnetic fields and time-scales all the ring diagrams are 
of the same order, and their resummation is necessary to 
obtain quantitatively correct results. 



and in the analogous way the summation over the odd 
orders in Eq. ( 68 1 gives 



N 

W(t) odd = JJexp (i[Xi{t) - arctanA^)]) , (73) 



and for FID we have Wp IU — W ( 



X Wodd- 

It remains to be shown that the calculation of W(t) 
using Eqs. (66 1 and (68 1 is feasible. The matrix with 



i and j indices labeling the nuclei has the dimension of 
N x N, which is too large for direct calculation. One 
possible solution is to use the continuum approximation, 
and replace the sums over the nuclei by integrals over the 
densities of hf couplings p(A) defined in Eq. (15 1. As for 



the errors, since we are effectively replacing the smallest 
possible Aij ~ A/N 2 with zero, we can expect these to 
arise at time t > N 2 /A ~ 1 s, which is an irrelevantly 
long time-scale. However, a drawback of this method 



is that the main merit of the approach from Eq. ( 66 1 



which is being able to obtain all the ring diagrams from 
a single calculation of the T-matrix, is lost. Only under 
certain conditions we will be able to obtain formulas for 
all the R k from the multiple integrals in the continuum 
approximation. 



We can take full advantage of Eq. ( 66 ) if we notice that 
Tki can be written as 



T u = fkiB{A k ,Ai) 



(74) 



where B(A k , Ai) is the coupling constant ( i.e. AiA 2 /4il 
for the lowest-order S^-conditioned interaction), and fki 
depends only on 



A w = (A k -A t )/2 

^>kl = ^k — . 



(75) 
(76) 



Now, starting from the continuum approximation for R k , 
we derive a formula just like Eq. ( 66 1 , but with an ef- 
fective T matrix replacing the original T. We coarse- 
grain the distribution p(A), dividing the relevant range 
of Ai into Ma slices, and deal with a T matrix of size 
M = MaNj, where Nj is the number of the nuclear 
species. We take a discrete set of A ak values for each 
species a with A ak = kA a with A a = A™ ax /M A . The 
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coarse-grained matrix is defined as 

T a k,pi = fki\n a nf3 I p a {M) j Pp{M) 



pi 



1/2 



B 2 (Ai,A 2 )dA 1 dA 2 



(77) 



Note that in the Ma — » oo limit Tki does not approach 
the original T-matrix, but the trace of the product of 
coarse-grained matrices approaches the continuum ap- 
proximation to the expression involving the original Ty. 
This definition of Tki is convenient because for Ma — 1 
it is equivalent to expanding the continuum approxima- 
tion expression to the lowest order in time (assuming 

While in the original T-matrix we had diagonal ele- 
ments equal to zero, in the coarse-grained matrix we have 
Tkki^O, because these matrix elements represents the in- 
teractions between ~ N /Ma nuclei, and removing them 
would correspond to an artificial lower bound on allowed 
Aij. Similarly as in the case of continuum approxima- 
tion, this approach is expected to lead to visible errors 
for t > NMa/A. It will turn out that the calculation 
of W{t) quickly converges at the time-scales of interest 
as we increase M (i.e. make the distribution of Ai more 
fine-grained), and it is enough to use M<C-/V. 

V. FREE INDUCTION DECAY 

We consider now the FID experiment, characterized by 
f(t; t) — 1 for r e [0, t] and otherwise. For the lowest 
order S^-conditioned interaction from Eq. ( 20 ) we have 



Tij = (1 - 5ij)y/a[a] 



AiAj 



40 A?- - ujfj 



icoki cos Aijt + Aij sin Aut 



The matrix elements for higher-order two-spin inter- 
actions (e.g. the ^-independent pair interaction from 
Eq. (22 1) are suppressed by powers of very small quantity 



A/NU 

When ujij ^> Aij (which is true in GaAs with N 



10 b and B~0.1T), the T i} from Eq. (78) for i and 



10 5 

j nuclei belonging to the same species are bigger than the 
hetero-nuclear matrix elements by a factor of the order of 
Wij/Aij. In this case, we can keep only the homo- nuclear 
terms in the T-matrix: 

A i Aj2smA ij t ^ 



Ti, 



(1 - %)(1 - 5 a [i].f}[j])a a [ 



A, 



This amounts to treating Nj nuclear systems as discon- 
nected, and to factoring the ring diagrams R 2 k into three 
contributions: Rk ~ ^2a=i -^fc ■ From this the factoriza- 
tion of W(t) follows: 



IF, 



FID 



e ''FID 



(*)> 



(80) 



with Wp^Ut) given by Eqs. (68), Rk calculated according 
to Eq. (166 b using the T-matrix from Eq. ( 79 1 . The clcc- 



(2) 

tron frequency shift A£l a comes from the second term 



in Eq. ( 19 1. Strictly speaking this term is a constant only 



for J = 1/2, but we will treat it in mean-field approxi- 
mation, i.e. we will replace it with its expectation value 
with respect to the narrowed state \J): 



Any 



A 2 

2Nn 



(81) 



Exact resummation for short times 



As a check of the calculation of W(t) using the matrices 
following from discretization of the distribution of the 
hf couplings, let us consider the limit in which we can 
calculate all the Rk analytically. We restrict ourselves to 
the short-times (A^t w At/N <C 1), which corresponds 
to t < 20 ps for a GaAs dot with N = 10 6 . Then we 
can approximate the homo-nuclear T-matrix for nuclear 
species a by 



rpo 



(1 - s ij)acx ji ^j Lt > 



(82) 



Using the continuum approximation we arrive at 

R k - f5!n a / A 2 Pa (A)dA) , (83) 



2fl 



which, using the density of hf couplings for a lateral dot 
from Eq. ( 16 ) gives us 



Rf.= 



2NQ 



(Vat? 



(84) 



Using the last equality we can resum the exponential se- 
ries in Eqs. (68 1 in the same manner in which we derived 



(78) Eqs. (721 and (73). Using Eq. (80), and noticing that 



Afl a — rj ai we arrive at the short-time expression for 
single-spin FID decoherence 



TF, 



FID 



-i arctan rj a t 



(85) 



The modulus of this expression was given in Ref. 1351 
Eq. (851, together with the factorization formula in 
Eq. ( 80 1 , is the main analytical result pertaining to the 



FID decoherence in III-V quantum dots for short times. 
Taking advantage of the long-range nature of the hf- 
mediated interaction, we have resummed the whole per- 
turbation series for W(t) obtaining a simple analytical 
expression. The characteristic FID decoherence time, de- 
fined by \W{T 2 )\ = l/e, is 



T 



2. short 



~A? 



(86) 



The same val ue was predicted within the PCA 
calculation ) 25 ! 26 ! but the form of the decay was different 
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there for t<^N/A. A similar characteristic time was a lso 
obtained using an equations of motion approach! 31 * 32 ! 

Let us stress that the T2 iS hort gives a characteristic time 
of decay only if T2 iSn ort is shorter than N/A, or equiv- 
alently when SI/ A < 1, i.e. for low magnetic fields. We 
comment on the issue of the high-field and long-time de- 
cay in the following Section. 



B. Relation to previous work on FID 



While our expression for W(t) is given by Eq. (68), the 
PCA result of Ref. [25] is given in our notation by 



Wfid,pca(£) = ex P 



1 



(87) 



The two formulas give indistinguishable results when 
R2 3> ^2^=3 Rk- I n t ne short-time regime this is fulfilled 
only when i?2^1- Then, both approaches give 
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FIG. 3: (Color online) Magnetic field dependence of the T2 
decay time characterizing the long-time exponential decay of 
coherence in FID experiment. The calculation is performed 
for Ino.5Gao.5As with iV = 10 4 and g e g = 0.5: the solid line is 
the T-matrix resummation of all the ring diagrams (M = 100 
was used), and the dashed line is the result of Eq. (94 1, with 
the total time given by T 2 = {J2 a T^ 1 ) -1 



which holds when rj a t <C 1, or equivalently t <C (N/A) ■ 
(SI /A). For B<1 T in GaAs we have SI/A<1, and this 
time-scale is shorter than the scale on which the Eq. ( 85 ) 



is valid. Thus, at low fields most of the coherence decay 



is well described by Eq. (85 1, with Gaussian decay being 



a good approximation only at very short times (when 
W(t) « 1). At low enough fields it should be possible to 
observe the regime in which rj a t 3> 1, and the coherence 
signal is given by 



www » n ■ 



Vat 



(89) 



which for GaAs means W(t) ~ l/t 3 . We stress that this 
formula holds when t > (SI /A) ■ (N/A) and t < N/A 
(and for moderate magnetic field, since uj a p^>Aij is also 
required) . At longer times and smaller B fields we have to 
evaluate numerically the full expression for W(t) instead 
of using Eq. fl85] ). 

Interestingly, at long times (t 3> N/A) a very differ- 
ent form of W(t) is obtained. The expressions for ring 
diagrams are of the form 



pa k _fc 



sin A\it sin A^t sin A^t 



dA^AJ.^A,) 



A 12 A 23 "' Afei ' 
and for the two-spin diagram we immediately get 



(90) 



R2 =«t / Pa{Ax)dAx / Pa (A 1 )dA 2 



A 2 A? 
1 -sine 2 Ai,-i 



ASl 2 



n * > 
(91) 



where sine x = sin x/x. As discussed in Ref. 26 in the 
t^>N/A limit, when sinc 2 Ayi — > jS(Aij), we obtain 



2S1 2 



P 2 a (A)A 4 dA 



2t 



^ 2, long 



(92) 



If R2^J2T=3 R k a * these long times, then the coherence 
decay is given by 



W(t » N/A)^Y[cxp(-t/T^ ong ) 



(93) 



showing that the influence of the bath at long times can 
be treated in a Markovian approximation. It is inter- 
esting to note that while the characteristic decay time 
T2 )S hort at low fields depended on the electron wave func- 
tion only through N defined in Eq. ([9]), the long-time de- 
cay constant I^iong does depencP^on the shape of W(r), 
specifically the distribution p(A) of hf couplings. For our 
wave-function from Eq. (12 1 and the corresponding p(A) 
from Eq. (15) we obtain 



2, long 



ZirNSl 2 
3 2 n n a n jA 5, 



57T 2 

192 



-V 

108 J 



(94) 



The same result for |W(i)| at long times has been ob- 
tained within a very different approach (GME using hf- 
mediated interaction from Eq. (20 1 ) in Ref. 33 However, 
it was shown there that the result from Eq. (94 1 is valid 
only for A/SI -c 1. Indeed, we also find that the full 
T-matrix calculation (i.e. resummation of all the ring di- 
agrams) gives a long-time exponential decay with T^iong 
time given by Eq. (94 1 only when A /SI -C 1. At lower 
fields we still htwe\W(t)\ ~ exp(-t/T^ ng ), but with 
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Taking larger than the value from Eq. (94 1. In Fig. pi we 



present the results of the calculation of magnetic fiela de- 
pendence of X^iong f° r Ino.5Gao.5As, where one can see 
that only for i?3>l T Eq. (94 1 gives a correct result. 



Let us stress, that the exponential decay dominates the 
coherence dynamics at high fields (f2 > .4) , at which the 
initial decay of | W(£)| in the short-time regime (discussed 
in the previous Section) is small, and most of decoher- 
ence occurs in the long-time regime. At low B most of 
the decay occurs in the short-time regime, and X^'dug 
shown in Fig. [3] at these fields is the characteristic time 
of exponential decay which emerges as long times when 
the coherence is already small. However, as we discuss in 
Sections |V C| and |VIII| we can safely argue that the ring 
diagram resummation approach is valid at low B fields 
only for short times. Thus, the low-B and long time re- 
sults, such as shown in Fig. [3] for B < 3 T, should be 
viewed with caution. 

To summarize, at low B our theory for FID predicts 
exponential decay at long times when |W(£)| is already 
very small, and the accuracy of this prediction has little 
practical bearing, while at high B our results agree with 
other approaches!^ For the purpose of quantum compu- 
tation using spin qubits, which is the main motivation for 
our work, the long-time decay is of no significance what- 
soever since one is only interested in the regime where 
quantum coherence is very high. However, for the pur- 
pose of establishing a closer connection between differ- 
ent theoretical approaches, the question of long-time and 
low-B behavior of FID is an important one, and we leave 
it for future investigation. Finally, let us mention that 
based on comparison of ring diagram theory with exact 
simulation of a small system we believe that for SE our 
theory works well also for long times and relatively small 
magnetic fields p3 



C. FID for uniform hf couplings 

We can test the ring diagram solution by applying it 
to the model system in which all the hf couplings Ai 
are the same (Ai — A/N), i.e. the electron wave-function 
is assu med to be constant inside of the dot and zero 
outside JMMII] We will also assume a homo- nuclear sys- 
tem with all the nuclei having the Zeeman splitting lo, 
and take the nuclear spin J— 1/2. In such a case, we can 
rewrite the second order term in the effective Hamilto- 



nian from Eq. ( 19 I as 



= ns z + Jj z + ^s z j z + s z ^^( ?:i ' 2 • 



:(J 2 -(J*) 2 ) ,(95) 



where 77 = A 2 /4Nn and J = uj - A 2 /4nN 2 , and J 
( J z ) are the operators of the total angular momentum 
of the nuclei (its projection on the z axis). It is natu- 
ral now to work in the basis of collective Dicke states 
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FIG. 4: (Color online) FID decoherence in a system with 
uniform hf couplings Ai. The solid lines are the real and 
imaginary parts of D m= o from Eq. (97 1, which is the exact 



solution for decoherence due to the lowest-order hf-mediated 
interaction. The symbols are results of calculation within the 
ring diagram theory. 



\Pjrn), in which j is the total spin of the nuclei, m is the 
eigenvalue of J z , and (3 is a permutation group quantum 
number PMl 

In order to obtain the narrowed state FID decay we 



plug in the above Hamiltonian into Eq. ( 29 ) in which the 



average corresponds to trace only over states with a fixed 
value of m. In the frame rotating with f2 frequency this 
gives us 



W m (t) 



exp 



x exp 



exp 



tt 2 

2n 



1 



N/2 

E 1 

j = \rn\ 



o 2 

n 

2n 



)t D m (t) 



(96) 



(97) 



where Q m = mA/N is the Overhauser shift, Z m — ^2 n j 
is the partition function in the narrowed state, and nj is 
the number of (3 states allowed for a given j (see Refs. [501 
I81I82I83II : 



N 
N/2-3 



N 

N/2-3-1 



(98) 



Note that Eq. ( 97 1 in fact predicts the coherence to revive 
at Poincare time tp = irN 2 fl/A 2 
where we used T2 short 



NT 2 . 

short short 3 

from Eq. (|S6 1 . Therefore this re- 



vival is inobservable in a realistic system, since other pro- 
cesses (such as spectral diffusion or even nuclear spin dif- 
fusion out of the dot) will alter the dynamics of the real 
system in the meantime and thus prevent the construc- 
tive interference of all the phases for different j states. 
On the other hand, the ring diagram solution in this 



1G 



case is given by 



W r \ nss (t) = e 



-iQ m texp[-iarctaiiT/t)] 



= „-i€l m t 



y/l + V 2 t 2 

^rings [J*) 



(99) 



Note that for uniform hf couplings the above solution is 
valid at all times (as long as we only consider the 
interaction), so that there is no transition to exponen- 
tial decay at long times. Furthermore, within the ring 
diagram approach we do not discern different m states 
(apart from the Overhauser shift f2 m ): in the diagram- 
matic derivation we have assumed that all typical states 
(with m<y/N) give the same W(t) in the frame rotating 
with f2 + Q m frequency. In such a frame the expression 
from Eq. ( f97| still has the exp(— itCl^/Q,) phase factor 
which is absent in the ring diagram expression. Thus, 
the shifts of precession frequencies of the electron spin 
are different in the two calculations. Practically, this is 
not an issue since both can be identified with the exper- 
imentally observed frequency. The question is whether 
the non-trivial parts of the decoherence function, i.e. D m 
and -D r j ngs agree with each other. Their comparison is 
shown in Fig. |4j where we plot D m= o calculated numeri- 
cally for iV = 2000 nuclei, showing a very good agreement 
between the two calculations. This results supports our 
claim that the ring diagram resummation is a good solu- 
tion for short times in realistic dots, since for t<^N/A 
the exact shape of the wavefunction does not matter, and 
all the diagrams depend only on TV defined in Eq. ^ . We 
have also checked that as long as m < VN, the shape of 
D m remains practically unchanged at time-scale shown in 
Fig. [4] confirming the result from Refs. I2H 26 that apart 
from different frequency shifts, all these narrowed states 
exhibit basically the same decoherence dynamics. 



D. Results for GaAs and InGaAs 

In GaAs we have three nuclear species, ^^Ga, 74 Ga, 
and 75 As, each having spin J — 3/2. The numbers of nu- 
clei per unit cell are given by n a — 0.604, 0.396, and 1, 
respectively, and the hf energies A a are given in Table [T] 
In Fig. [5^, we show |Wp ID (i)| for various N and B calcu- 
lated using our theory, and in Fig. [5 b-c we present their 
comparison with the PCA calculation!^!^ In Fig.[5]3 one 
can see that at low B, when the decay occurs at times 
shorter than N/A, there are visible differences between 
the fully resummed solution and the use of only two-spin 
ring diagram. At higher B used in Fig. pj;, when the 
decay occurs at longer times, our calculation is in much 
closer agreement with the PCA results of Ref. [25j since 
in the high B and long-time regime the only difference 
between the two results is the presence of corrections to 
T2,iong in our theory (which are quite small at B = 10 
T in GaAs). The additional phase dynamics due to the 
ring diagram resummation is shown in Fig. [6] where real 
and imaginary parts of Wp ID are plotted. 
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FIG. 5: (Color online) Narrowed state \W s {t)\ for FID cal- 
culated for two GaAs dot of different sizes and various mag- 
netic fields, (a) Our theory calculated with coarse-grained 
T-matrix with Ma = 10. (a) Red (blue) lines correspond 
to N = 10 5 (10 6 ), while solid (dashed) lin es correspond to 
5 = 0.1(1) T (g ft = — 0.44 is used). Eq. (85 I corresponding to 
Ma ~ 1 gives practically the same results at this time-scale, 
(b) Comparison of our result with all the ring diagrams re- 
summed (solid line) with the pair correlation approximation 
(PCA) solution (circles) for N = 1Q 6 and 5 = 0.1 T. The be- 
havior of I W(f)| after the initial decay (for t> 1 /is) is ~ l/t 3 
and ~ exp(— t 2 ), respectively, (c) Analogous comparison for 
TV = 10 5 and B — 10 T, when the decay occurs in the long 
time regime (t^$> N/A), and the pair approximation is much 
closer to our result. 



The optical experiments on an ensemble of In- 
GaAs dots have been intcrpretecPSl as measurement of 
narrowed-state FID (note, however, that the physics of 
the nuclei influenced by a train of optical pulses is richer 
and is still being investigated, see Refs. 154184"]) . The ob- 
served characteristic decay time was between 3 /is for 
larger,^ and 0.2 /is for smaller dots. 84 In Fig. [7] we 
present our results for FID decay in In^Gai^^As dot with 
N = 10 4 — 10 5 nuclei and indium concentration of a; = 0.5. 
The parameters of the Hamiltonian are given in Table [T] 
the concentrations of 113 In and 115 In are n a — 0.0428a; 
and 0.9572a;, respectively, and both of these nuclei have 
spins J = 9/2. Note that for N = 10 4 practically the 
whole decay occurs in the long-time Markovian regime of 
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FIG. 6: (Color online) Comparison between the numerical T- 
matrix calculation of Wp ID (i) and the analytical approxima- 
tion from Eq. ( |85[ ). The solid lines are the real and imaginary 
part of W s (t), and the symbols are from Eq. (851. The cal- 
culation is done for GaAs with N = 10 5 at B = 1 T. In the 
coarse-grained T-matrix calculation we have used Ma =30. 
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FIG. 7: (Color online) FID decay in Ino.5Gao.5As for A = 10 4 
(blue) and 10 5 (red), at B = 5 T and with g e ff = 0.5. The solid 
lines are | W(£)| , while the dashed (dot-dashed) lines are the 
S x {S y ) components of the spin normalized to unity. We have 
assumed that at t — the electron spin is initialized along the 
x axis. 



exponential decay, but at the field of 5 T we expect the 
long-time solution to be reliable The decay is practi- 
cally completely dominated by interaction with In nuclei: 
due to their spin of 9 /2 they are much more efficient at 
decohereing the electron than ,7 = 3/2 nuclei of Ga and 
As. We also remark that using the results from Ref. UHl 
we have estimated the spectral diffusion decoherence time 
in such small InGaAs dots to be about 5-10 fis, justify- 
ing our concentration on decoherence due to hf-mediated 
interactions which occurs on sub- microsecond time-scale. 



VI. SPIN ECHO 

It was first observed in numerical simulations^ that 
at high B practically all of the decoherence due to hf- 
mediated interaction was removed by the SE sequence. 
This result was explained in intuitive and tra nspare nt 
way using the pseudospin approach in the PCAPSEH In 
these articles, the lowest-order S^-conditioned interac- 
tion was treated together with the diagonal dipolar inter- 
action (—2bijJ?Jj in Eq. ([Hi), resulting in a near-perfect 
reversal of the nuclear spin aynamics by the SE sequence. 
The magnetic-field dependence of the small remaining de- 
coherence was very wealP^ in the high B regime consid- 
ered there. It is crucial to note that in these papers the 
nuclear bath was assumed to be either homonuclearj^or 
the interactions betw een n uclei of different species were 
(justifiably) neglecteoPMl at BwlO T. 

The effect of the SE on decoherence due to the dipo- 
lar interactions is much less dramatic, and it is now 
well established that in the very high magnetic field 
limit th e SE decoh erence is purely due to spectral 
diffusion ! 25 ! 26 ! 27 ! 2 -^^ Here we consider much lower mag- 
netic fields and concentrate only on the role of the hf- 
mediated interactions. 

The following observation offers the key insight into 
which interaction channel is important in this case. If 
we remove all the nuclear intrabath interactions with ex- 



ception of the S z -conditioned terms (i.e. V\ in Eq. (23)), 
and treat these remaining terms in the secular approx- 
imation (allowing only for processes conserving the nu- 
clear Zeeman energy), then there is no SE decoherence, 
i.e. Wse(£) = 1- Indeed, the nuclear interaction Hamilto- 
nian in the secular approximation commutes with nuclear 
Zeeman energy Hz n , and we have 



-m±t 



-i(H Zrl ±H A ±V 1 )t 



e -iH Zn t e ±i{H A + V x )t 



(100) 



so that the product of evolution operators from Eq. ( 30 1 



gives unity in Eq. ( 26 1 . The same holds for any balanced 
dynamical decoupling sequence. Thus, there are three 
ways of getting finite SE decoherence from hf-mediated 
interactions. First is to reintroduce the dipolar inter- 
actions and consider the terms in the perturbation the- 
ory which mix them with the hf-mediated interaction (as 
in Refs. H5J.26J. Second is to consider the non-secular 
processes, and the the third way is to consider the S z - 
independent hf-mediated interactions. 

Here we consider the two latter approaches for the fol- 
lowing reasons. Mixing of the hf-mediated and dipolar 
interactions creates terms in the effective Hamiltonian 
which are partially local, and the number of terms in ex- 
pansion of W(t) due to these interactions does not scale 
as N k with the number k of spins involved. Furthermore, 
a ring diagram (having all the nuclei distinct) with at 
least one S^-conditioned interaction will still be equal to 
zero for SE. In order for the term in perturbation ex- 
pansion to be non-zero, one has to consider the cases 
of repeating indices, e.g. by doing the full diagrammatic 
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calculation as in Ref. [3D] The need to repeat the nuclear 
index in the summation brings us back to the "term- 
counting" argument. Finally, separating the dipolar and 
hf-mediated interactions allows for intuitive and trans- 
parent treatment: we can deal with the dipolar interac- 
tion using the well-tested cluster expansion theory, ^ 28 * 29 ! 
and treat the hf-mediated interactions separately using 
the ring diagram resummation which is appropriate for 
the case of long-range interactions. 



A. T-matrix solution for Spin Echo 



We repeat the calculations from Sec. IV using f(t; r) 
for the SE is shown in Fig. [T]d. For the two-spin S z - 
conditioned interaction from Eq. (|20[) we get 



8zw fe/ e^ fc! * /2 sin Afc 't^' tsin Afc '7^' t 



A 2 

A kl 



• CO 



kl 



(101) 

As expected from the preceding discussion, we have 
Tki = when the two nuclei belong to the same species 
(i.e. when ujki — 0). At moderate B fields (when uj k i 3> 
Aki) and at short times we obtain for hetero- nuclear pairs 



-e ZUJklt/2 sin 2 



(102) 



The coarse-grained T-matrix from Eq. ( 77 1 with Mj = 1 
(expected to lead to accurate results at t 
is then given by 



T a/ 3 = (1 - 5af3)^/a a af3^/n a n0 
- 21 e-^sin 2 ^ 



lese short times) 
(103) 



With this we can write the ring diagram contribution as 



Nj 



R 2k (t) 



A 



(104) 



where A; are the eigenvalues of the T matrix. As in 
Eq. (72) we arrive then at the closed expression for the 



dccohcrcncc function due to the non-secular flip-flop pro- 
cesses: 



N, 



IF 



SE 



non-sec 



i 



(105) 



=i Jl + A 2 



The case of Nj = 3 is experimentally relevant for GaAs 
quantum dots. There we obtain A^(i) = 0, ±A(i) which 
leads to 



1 



(106) 



where we have identified R 2 = 2X 2 , and A 2 is given by 



A 2 (i) = f 12 {t) + f 13 (t) + f 23 (t) 
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FIG. 8: (Color online) Spin echo decoherence WsE(i) in GaAs 
at low B fields. The parameters of GaAs are the same as used 
in FID calculations. Full T-matrix calculation was performed, 
but the results for N = 10 6 (solid lines) are indistinguishable 
from the approximation of Eq. (106 1, while some differences 



between this analytical approximation and the numerics arise 
for the smaller dot with N — 10° (dashed line) at longer times 
(see Ref. The results for 25 mT and 50 mT show the 

saturation of decay at fields smaller than B c (see the text). 
The comparison of solid and dashed red lines (both for 0.1 
T, but different N) illustrates how the \ow-B decay scales 
with the dot size. The dots are the PCA result for 25 mT 
and N = 10 6 , which should be compared with the solid blue 
line showing the result of calculation of all the ring diagrams. 
Similarly as in case of FID in Fig.|5]D, PCA overestimates the 
decay beyond its initial stage. 



and T a p are given in Eq. ( 102 1 



Let us mention that analogously to the FID case, the 
higher-order two-spin interactions are giving negligible 
contributions compared to the one discussed above. For 
example, when we use the S' 2 -indcpendcnt term from 
Eq. (22 1 as the interaction, we get the formula analo- 
gous to Eq. (101 1, but with the multiplicative factors 
of ujki in the numerator replaced by Aki and with Bki 
replaced by C a = A k Ai(A k + Ai)/I6fl 2 . Now both sec- 
ular and non-secular terms are non-zero, but the secular 
terms are much larger for uj k i 3> A k i, and the calcula- 
tion of this contribution to decoherence parallels the one 
performed for FID in Sec. |Vj The result is that the cor- 
rections due to this interaction are completely negligible 
compared to the decoherence due to the lowest-order in- 
teraction for magnetic fields larger than 10 mT in a dot 
with 7V=10 5 - 10 6 nuclei. 



B. Spin Echo results and discussion 

At moderate B fields the SE decoherence is fully de- 
scribed by the two-spin ring diagram R 2 appearing as a 
result of resummation in the denominator of Eq. (1061. 



(107) For GaAs R 2 (t) is given by a sum of three positive pe 
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riodic functions, each having an amplitude ~ B~ A . Let 
us define the magnetic field B c for which these ampli- 
tudes become of the order of 1 . This B c corresponds to 
the electron Zecman splitting of fl c ~ ^/tA/VN, where 
r = Q/u), with lu being the typical difference of the nu- 
clear Zeeman splittings, is of the order of 10 3 . Then, for 
B > B c Eq. ( 106 ) adds only a small oscillatory compo- 
nent to the total W(t), and the actual SE decay is due 
to the spectral diffusion, which occurs on time-scale of 
more than 10 fjs for dots with N> 10 5 , with the actual 
value depending on the shape of the dot J 28 * 29 ! 

On the other hand, at B < B c (corresponding to B < 
250(75) mT in GaAs with | 5cff | = 0.44 and N= 10 5 (10 6 )), 
the decay due to the lowest-order hf-mediated interac- 
tion is significant. J R 2 (i) is a sum of periodic functions, 
but their respective periods are incommensurate, and it 
is not expected for R%(t) to come back to (for W(t) to 
come back to 1) at a finite time t. Only partial "rephas- 
ing" occurs, and at B < B c the minimal value of i?2(t) 
becomes large and the time at which it is first achieved 
becomes long, so that the predicted SE decay is then from 
practical point of view irreversible, since spectral diffu- 
sion makes W{t) decay to zero at a time-scale of > 10 //s 
anyway. The characteristic time of SE signal decay due 
to the hf-mediated interactions for B <B C is given by 



SE 



N 



3VFT* , 



(108) 



which in GaAs translates to Tse~ IOOT2 . These results 
arc illustrated in Fig. |] for GaAs dots with N= 10 6 and 
10 5 nuclei. Eq. ( 108| predicts the characteristic low-B 



(10 fa ) 



decay time of « 0.5 fis (1.5 /its) for N = 10 b 
agreement with the results shown in Fig. [8] 

The calculated decay at low B (when the signals cease 
to depend on the magnetic field) is in qualitative agree- 
ment with Tse ~ 0.3 — 0.4 /is measured in Ref. [12] for 
B « 50 — 70 mT. The prediction of our theory is that at 
B ten times larger the decay will be incomplete, having 
oscillations on a microsecond time-scale with a period 
proportional to B. Let us note that according to the 
theory presented here, the increase of Tgg with increas- 
ing B becomes visible only once B becomes larger than 
B c = il c /g c ff/^B, and it should be accompanied by an in- 
complete decay at longer times. At even higher B the 
oscillations of the SE should become visible. Let us note 
that it might be possible that the decay seen in Ref. H2l 
was in fact incomplete (like the line line corresponding 
to 5 = 0.1 T and N= 10 6 in Fig. [§), since the measure- 
ments were done in a time window not much longer that 
the observed Tse, and the exact value of signal (i.e. cur- 
rent through the dot) corresponding to zero coherence 
could be uncertain. 

In Fig. [9] we show a comparison of W(t) for the single- 
spin FID and SE calculated for the same GaAs dot. As 
expected,^ the coherence time is increased in the SE ex- 
periment compared to the single- spin FID. Note however 
that this is not entirely trivial, since the SE decay is dom- 
inated by processes of different character than the ones 
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FIG. 9: (Color online) Comparison between the single-spin 
FID decay and SE decoherence in GaAs quantum dot with 
7V = 10 5 and iV = 10 6 at B = 0.1 T. 



contributing to the FID dynamics. 



VII. RESULTS FOR DYNAMICAL 
DECOUPLING WITH MULTIPLE tt PULSES 

The calculation from the previous Section can 
be easil y genera lized to dynamical decoupling (DD) 
sequence d 57 1 58159 1 of many ideal 7r pulses driving the elec- 
tron spin. DD with various pulse sequences was con- 
sidered theoretically for a spin bath both in a high B 
and long ti me regime wh en the dipolar interactions play 
a large rolej 26 * 60 * 61 * 62 * 63 ! and f or ve ry low B fields when 
one has to use exact numerics! 41 * 42 ! Here we focus on the 
heteronuclear system at low B considered in Sec. |VI| 
and calculate the W(t) for an experimentally realistic 
case of applying a few pulses. As in Sec. |VI[ we use a 
short-time approximation, so that only few first microsec- 
onds of coherence evolution are considered for GaAs with 
N ~ 10 6 nuclei, and we assume a moderate B field at 
which Lu a /3^>A/N. Only lowest-order hf-mediated inter- 
action is considered, since it is the dominant source of 
the SE decay in this situation. 



At short times we have W(t) given by Eqs. ( 106 )-( 107 1, 
only a different filter function f(t; r) has to be used in 
calculation of T Q| g. It is convenient to rewrite R2(t) as 



R 2 = c ^ f dT i I dT *f& n)e^ (Tl " T2) /(*; t 2 ) 
„ «j„ Jo Jo 

(109) 
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Let us introduce the "spectral density" function 
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FIG. 10: (Color online) Decoherence under various dynam- 
ical decoupling pulse sequences in GaAs quantum dot with 
N = 10 6 at 6 = 0.25 T. The results for SE, 2-pulse CPMG 
sequence, and 5-pulse CPMG, CDD, and UDD sequences are 
shown. 



and also the filter function in the frequency domain de- 
fined by 



F(ujt) 



2 _ 2 



(112) 



where f(t;u>) is the Fourier transform of f(t; r) with re- 
spect to t. Then we can rewrite Eq. ( 109 1 as 



Ri = 2 



2^ M ^ 



(113) 



The purpose of this formulation is the following. The 
pure dephasing decoherence due to S z /3(t) coupling to 
classical Gaussian noise (3{t) having spectral density S(u>) 
is gi ven by W(t) — exp[— R 2 (t)/2] with R 2 given by 
Eq. \IVZ). For derivation see e.g. Refs. 1711721 Thus, at 



very short times, when R% <C 1, the decoherence due to hf- 
mediated interactions between different nuclear species in 
regime of interest can be mapped onto decoherence due 
to classical Gaussian noise. It is interesting to note that 
a similar formal equivalence between initial decoherence 
of a quantum model under pulses sequences and classical 
noise problem has been derived for a very different bath 
Hamiltonian in Ref. [701 Furthermore, for the case at 



hand we have Eq. ( 106 1 , which allows us to use Eq. (113) 



to describe the decoherence at longer times, at which the 
analogy to the classical Gaussian noise does not hold. 

The frequency domain filter functions F(ujt) were 
calculated in Ref. [72] for different pulse sequences 
aimed at preventing pure dephasing, namely the CPMG 
sequence,^ the concatenations of spin echoJSMLLJ^ 
the peri odic seq uence of it pulses,^ and the Uhrig's 
sequence! 59 1 63 1 85 1 For example, for SE we have F(z) — 
8sin 4 [z/4], which leads us back to the results obtained in 
Sec. |VI| The two pulses CPMG, concatenated spin echo 



(CDD), and Uhrig's sequence (UDD) all correspond to 
t/4 — 7T — t/2 — 7T — i/4 sequence (see Fig. [T|d), with the 
filter 



F(z) = 128 cos 2 - sin 6 



(114) 



and F(z) functions for sequences with more pulses are 
given in Ref. [751 

The common feature of all the DD sequences is that 
their F(ujt) functions are filtering out large part of S(oj) 
in Eq. (113 1 at frequencies smaller than where n 

is the number of pulses. In fact, the Uhrig's sequence 
(UDD) is most efficient at this task when th e noise spec- 
trum has an upper frequency cutoff ) 63 ^ 72 ^ 86 ^ which is ex- 
actly the case here (the cutoff being the maximal |w a /3|)- 

In Figure 10 we present the results obtained for W(t) 
in GaAs for DD sequences with up to 5 pulses. As ex- 
pected, the UDD sequence makes W(t) very close to 1 
for the longest time. The characteristic features visible 
for times at which (W(t) is smaller, i.e. the peak in signal 
for five pulse CPMG sequence and a plateau for the CDD 
sequence, originate from the discrete nature of the noise 
spectrum. The peak in W(t) for n = 5 CPMG comes 
from the io a pt factor (with u> a ^ frequency for 71 Ga and 
75 As nuclei) matching a minimum of the F(wt) filter at 
uit~ 37r. The plateau in the CDD signal has similar ori- 
gin. In fact, such features have been observed in a recent 
study of dynamical decoupling from classical noise in ion 
trap qubits, 87 where a plateau in W(t) was related to the 
presence of sharp peak in S(uj). 



VIII. CORRECTIONS DUE TO HIGHER 
ORDER TERMS IN THE EFFECTIVE 
HAMILTONIAN 

In the preceding sections we have presented calcula- 
tions of decoherence due to hf-mediated interactions in- 
volving pairs of spins. Among such interactions in the 
effective Hamiltonian H, the most important one was 
the lowest-order S^-conditioned one from Eq. (20 1, and 



the contribution of higher-order two-spin interactions was 
shown to be very small using the ring diagram approach. 
What remains to be shown is that the influence of higher- 
order multi-spin interactions (such as the three-spin one 
in Eq. (21 1) is small at the time-scale of coherence decay. 

First, let us make a general remark about the con- 
vergence of the operator series H^ n \ If one applies a 
commonly used operator norm (see e.g. Ref. [ 



max I ($1 A 1$) 

($|*)=i 



(115) 



to the leading H^ n ' terms (i.e. terms containing n spin 
operators), one gets an estimate 



\Jj{n)\ 
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(116) 
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so that it seems that the expansion is convergent only 
when ^4/0 < 1. However, the state |$) which maximizes 
Eq. (1151 is a highly entangled state of nuclear spins. 



Here we are interested in evolution of nuclear spin sys- 
tem averaged using an uncorrelated thermal ensemble, 
and such highly entangled states do not appear in the 
calculations. Thus, it is plausible that the condition for 
convergence (in the sense of giving us well defined results 
for W(t)) of the expansion of H, is much less restrictive 
than A/CK1. 

The full treatment of decoherence due to leading order 
terms in H would require a derivation of a diagrammatic 
technique involving multi-spin interactions. Instead of 
pursuing this path, we choose to give a more limited, but 
practically sufficient argument. We calculate the lowest 
order contribution to W(t) due to a multi-spin interac- 
tion, e.g. for three-spin (3s) interaction from Eq. (21 ) we 
obtain 



(2) 
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and although the analogous expressions for due to 

interactions involving more than three spins get more 
complicated, they are suppressed by higher powers of S 

(2) 

compared to ■ In fact, on the time-scale of t<^N/A, 

(2) 

while W^ s can be larger than 1 (thus forcing us to resum 
the higher-order terms), for <5<C 1 the lowest-order contri- 
butions from multi-spin interactions are still negligible, 
and it is not necessary to go to the higher orders of per- 
turbation theory to conclude that these interactions are 
irrelevant on this time-scale. Note that the key results of 
this paper were obtained for magnetic fields at which the 
coherence was decaying significantly in the short time 
regime. Thus, we arrive at conclusion that as long as 
5 <C 1 , the contribution of multi-spin interactions to the 
full W(t) is very small on the time-scale on which the 
coherence decay due to the two-spin interaction occurs at 
low fields. The question of convergence of our approach 
at longer times, relevant for smaller N or larger B (but 
not so large as to give A/£1<1, when the convergence of 
the effective Hamiltonian expansion becomes more obvi- 
ous), is left for future investigation. 

IX. CONCLUSION 

We have presented a detailed description of the theory 
of pure dephasing decoherence of an electron spin cou- 
pled to a nuclear bath by hyperfine interaction. At finite 
magnetic field B one can perform a canonical transforma- 
tion which eliminates the direct electron-nuclear spin-flip 
terms in hf interaction in favor of an effective pure de- 
phasing Hamiltonian containing hf-mediated long-range 
interactions between the nuclei. We have argued that 
when the bath is thermal (uncorrelated and unpolarized) , 
one can use the lowest order hf-mediated interaction to 
calculate the spin decoherence as long as the electron 



spin splitting is larger than the rms of the Overhauser 
field distribution ~ A/VN (this correspond to B > 3 
mT in GaAs dots with 10 6 nuclei). The solution is pos- 
sible due to long-range nature of the hf-mediated inter- 
actions, which couple with comparable strength all the 
N nuclei within the bulk of the electron's wave-function 
(N being larger than 10 4 in III-V based quantum dots). 
In such a case we can identify the leading (in terms of 
1/N expansion) diagrams in the perturbative series for 
the decoherence function, and then resum all these ring 
diagrams. At short times (when t -C N/A, where A is 
the hf coupling energy) we can perform this resummation 
analytically and obtain closed formulas for decoherence 
in Free Induction Decay (FID), Spin Echo (SE), or any 
other dynamical decoupling sequence of pulses driving 
the electron spin. At longer times, the solution can be 
obtained numerically at very low computational cost. 

Compared to previo us work o n decoherence due to hf- 
mediated interactions l 25 l 26 l 3 ° l 33 l we have accounted for 
all the leading terms in the diagrammatic expansion of 
decoherence time-evolution function. This allows us to 
firmly establish the regimes in which the pair correlation 
approximation^^!] is quantitatively applicable, and also 
to make comparisons with results obtained using a very 
different Generalized Master Equation approach The 
resummation of all the ring diagrams leads to modified 
results for single-spin FID (i.e. FID measured in a nar- 
rowed nuclear state with well-defined polarization) decay 
at low magnetic fields comp ared t o the pair-correlation 
approximation (PC A) results w hile at high fields we 
recover the PC A solution given in these papers. We also 
found nontrivial phase dynamics of the spin coherence 
in the FID experiment. The ring-diagram approach to 
FID was tested by showing that it agrees with an exact 
analytical solution in the case of uniform hf couplings. 

For SE decoherence, apart from including all the ring 
diagrams in the calculation, we have identified the most 
important mechanism of SE decay in III-V quantum dots 
at low magnetic fields: the hf-mediated interaction be- 
tween spins of different nuclear species. While at high 
B this interaction channel leads only to a small mod- 
ulation of the SE signal, at low B used in recent SE 
experimentation GaAs dots this leads to coherence decay 
on a sub- microsecond time-scale. This result is in quali- 
tative agreement with the measurement from Ref.[T2] We 
have also given predictions for the effect of dynamical de- 
coupling pulse sequences on decoherence in the regime in 
which the SE has been measured. 

We believe that this work extends the analytical theory 
of spin decoherence due to interaction with the nuclei 
down to magnetic fields much smaller than previously 
assumed, and leaves only the highly nontrivial zero field 
limit to exact numerical approaches. Using the theory 
presented here one can make qualitative and quantitative 
predictions, particularly in the experimentally relevant 
regimes of low magnetic fields and relatively short time 
scales. 
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Now we expand iS into operators of increasing order in 



APPENDIX A: CANONICAL 
TRANSFORMATION LEADING TO THE 
EFFECTIVE PURE-DEPHASING 
HAMILTONIAN 

Following Rcf. 88 we divide the Hilbcrt space of our 
problem into two parts, A and B, spanned by eigenstates 
of H%, and we assume that for all eigenstates m e A 
and I £ B we have E^ — E?k&L, with f2 being much 
larger than the spread of energies within A and B and 
the energy scale of perturbations which we will consider. 

The full Hamiltonian is then written as 



H = Hn + Hi 



Ho 



Hn + H' 



(Al) 



where H — flS z , Hi has only matrix elements within 
the A and B subsets of states (i.e. the diagonal hf cou- 
pling AiJfS z , nuclear Zeeman, and dipolar interactions 
are contained in Hi), and H2 has only matrix elements 
between the states from blocks A and B (and actually 
H-2 = V s f, the spin-flip part of the hf interaction). We 
want to devise a transformation which lead to a new 
Hamiltonian H which has no matrix elements between 
the states from A and B subspaces. We write this Hamil- 
tonian as 



H 



* ' 7)1 



(A2) 



where S is an anti-Hcrmitian operator, and we have de- 
fined the nested commutator: 



[H,S]W = [...[[H,S},S], S] 



(A3) 



with S appearing n times. 

We refer to the operators having nonzero matrix ele- 
ments only within the A and B subspaces as block diag- 
onal (BD), and to those having nonzero matrix elements 
only between the states from A and B subspaces as block 
off-diagonal (BOD). A product of two BD or two BOD 
operators is a BD operator, while a product of a BD and 
a BOD operator is a BOD operator. In order to obtain 
a closed hierarchy of equations for S we need to assume 
that S is a BOD operator. Then the condition for the 
BOD part of H to vanish is 



00 00 

E (2^fTj!l*° + *'^l' 2 " +1) + E (i^ 1 * 2 ^ 1 ' 2 "' = 

n=0 v ; n=0 y ' 

(A4) 



^ (2n)\ 

n=0 y ' 



S =<?« +£(2) 



(A6) 



and by plugging this expansion into Eq. (A4) we obtain 
an infinite set of equations: 



[H ,S^] = -H 2 , 
[H ,sW] = -[Hi,sW] , 

[H ,S^] = -[HiJM] - ^[[H^S^JV] ,etc. 



(A7) 
(A8) 

(A9) 



Now we define the operator A: 

A = P A - P B 



(A10) 



in which Pa(b) is the projection on the subspace A(B). 
With the energy separation between the subspaces: 



n 



E° A 



E B 1 



(All) 



where E A B are the typical zeroth-order energies (i.e. the 
diagonal elements of H\ ) , we arrive at 



[H ,X] ps — fLAA 



(A12) 



for X being a purely block-off-diagonal operator (such as 
S). This approximation corresponds to neglecting terms 
of the order of A/vN and tOi compared to ft. This will 
affect the values of coupling constants in certain terms 
in H^ with n>2, but it will not alter the structure of 
multi-spin nuclear interactions in H and the scaling of 
their coupling energies with A and O. The formulas for 
the second-order effective Hamiltonian derived without 
making this approximation are given in e.g. Refs. 1331381 
Using the following properties of A: 

A 2 = l, [Hi,A}=0 , {H 2 ,A} = 0, (A13) 

we obtain the explicit formulas for S^: 

= ^H 2 A , (A14) 

sW = ±[Hi,H 2 ], (A15) 

^ (3) = ^[HidHuH 2 ]]A+^[[H 2 ,H 2 A],H 2 A]A 

(A16) 
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Using these, from Eq. (A5l we can H up to the fourth 



order in H' 






= 


Hi , 


(A17) 


H (2) = 


\[H 2 ,SW} = 


±[H 2 ,H 2 A] = ±H 2 2 A,(A18) 


H {3) = 




^[H 2 ,[H U H 2 }] , (A19) 


H^ = 


A 




= 




^ 2^ 2 Hi — H 2 HiH 2 H\ 




—HiH 2 HiH 2 


+ H 2 H\H 2 - HfjA . (A20) 


n hf 


- ^-4ft2 

i 


fz | \ " A l Aj(A l + Aj) p p 

1,3 



In higher orders the formulas get even more complicated. 
However, for n odd we always have expressions involving 
products of even number of H 2 and odd number of Hi , 
whereas for n even we get ffW being a sum of products 
of even numbers of Hi and H 2 operators, multiplied then 
by A. 



II B 



where 



The full expression for H^ is given in Sec. 
also certain terms from H^ were shown. Here we give 
all the terms in H^ which arise when we neglect the 
dipolar interactions in Hi (i.e. we keep only the Zeeman 
and diagonal hf term in Hi): 



-A 

I 

AH, 



AjAj(Ai + Aj) j + j 



16ft 2 



a 3 



EA 2 A 



i^3 



A l A j A k ~ + - ~ z 
1 * k 



i±3 



4ft 2 



(J-J+ - J+J-) 



(A21) 



The last term above is the three-spin interaction from 
Eq. ( 22 1 . Terms such as this one, which contain n nuclei 



in n-th order of expansion of H, axe potentially the most 
"dangerous" ones when one considers the convergence of 
the expansion. For example, in the the next order we 
have such a "maximally nonlocal" term given by 



h[ a ) — —S 



, AjAjAkAi - + - - + - 

2-~> J * J 3 J k J l 



(A22) 



APPENDIX B: REPLICA TRICK PROOF OF 
THE LINKED-CLUSTER (CUMULANT) 
EXPANSION 

The simplest proof of the linked-cluster theorem fol- 
lows from the replica trick.^1 The starting point is the 
expression for W(t) from Eq. (40 1. We then introduce 



a set of R replicas of the nuclear system (i.e. R inde- 
pendent quantum dots). Since the generalized (contour 
ordered) exponent still has the defining property that 
for independent (commuting) operators X 



-,X+Y . 



and Y, the R-ih power of the original W is simply equal 
to the W(t) calculated in the space of all the replicas. 
The Hamiltonian in this space is then simply the sum of 
the Hamiltonians corresponding to R uncoupled systems, 
and the interaction in Eq. (40 1 becomes V=Yl r r' ^r,r'V r , 
i.e. it is diagonal in replica indices r,r'. Now we use the 
fact that 



w R = e RlnW = 1 + RlnW+ y (R^W) m (B1) 

^ ml y ' 

m=2 



so that the part of W R which is linearly proportional to 
R is simply In W. When we look at the ring diagrams 
derived using the replicated interaction, we see that since 
the different replicas are independent, each linked (ring) 
diagram has to be summed over an additional label r. 
Due to the equivalence of different copies of the system, 
each ring diagram acquires a factor of R. Consequently, 
the diagrams which are products of k rings acquire a 
pref actor of R k . Therefore, In W is given by a sum of all 
the linked terms in the perturbation expansion of W. 
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